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Abstract 

We consider the short-distance behaviour of the product of the Noether 0{N) 
currents in the lattice nonlinear fi-model. We compare the numerical results with the 
predictions of the operator product expansion, using one-loop perturbative renorma- 
lization-group improved Wilson coefficients. We find that, even on quite small lattices 
(ma ~ 1/6), the perturbative operator product expansion describes that data with 
an error of 5-10% in a large window 2a ^ x ^ m~^. We present a detailed discussion 
of the possible systematic errors. 



1 Introduction 



The lattice regularization is nowadays the best theoretical tool for the study of nonper- 
turbative phenomena in field theories. It has been extensively applied to QCD, providing 
quantitative predictions for many non-perturbative quantities, and it has also been exten- 
sively used in the study of strong-interaction effects in weak processes. One of the basic 
problems one has to deal with in these investigations is the determination of renormalized 
operators. Referring to the simple case of a purely multiplicative renormalization, given 
a "bare" operator Olat on a lattice of size L, we want to find a renormalized operator 
related to Olat by 

ORifx) = Zo(A//i,Aa,aL-i) Olat- (l-l) 

In Eq. ( |1.1| ) we have emphasized the dependence of the renormalization constant Zo upon 
the various scales of the problem: the lattice spacing a, which encodes the dependence 
on the bare coupling constant Qq, the renormalization scale n, the physical scale A (the 
so called "A-parameter" ) , which breaks the scale invariance of the continuum theory, and 
the size of the lattice L. The constant Zq should be defined so that Or^h) has a finite 
continuum limit. In other words, the limit Aa — >■ (i.e. gi oo), L/a oo, at fixed A//i, 

0r(A//i)= lim^^lim (1.2) 

should be finite.[] It is important to notice that the two limits in Eq. ( |1.2|) do not commute, 
the infinite-volume limit should be taken before the continuum extrapolation Aa — 0. The 
renormalization constant admits a smooth infinite-volume limit at fixed a. The behaviour 
in the continuum limit is instead singular and it is predicted by the renormalization group 
(RG). An exception is represented by operators that are conserved, at least in the continuum 
limit. In this case, the renormalization constant has a smooth continuum limit and does not 
depend on the scale /x. Examples of such operators are the vector and axial flavour currents 
in QCD and the energy-momentum tensor. In this case, the renormalization constants can 
be found by requiring the validity of the Ward identities related to the conservation of the 
operator 

More difficult is the computation of the renormalization constants when the operator 
is not conserved. In this case, one can use lattice perturbation theory. However, these 
calculations are very tedious and long and one is usually restricted to one-loop order. A 
general non-perturbative strategy has been proposed in Ref. [^j and extensively studied 
since then, see, for instance, Ref. 0. It provides a direct method for the computation 
of the renormalization constant appearing in The main difficulty is the continuum 



limit, Eq. (L2), which requires to keep under control the lattice artifacts and the finite-size 



effects. In practice, the method works only if there exists a window 

L/a> — >1, L/a>-^>l, (1.3) 

jj,a Aa 

in which lattice artifacts 0{fia,Aa) and finite-size effects 0{a/L) are sufficiently small. 
Conditions (|1.3|) are not enough since we also want to relate (9™°* to renormalized operators 



^Eq. (^^) is a symbolic notation. Explicitly it indicates that correlations of Or{^) with renormalized 
fields have a finite continuum limit. 
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defined in continuum, say MS, perturbation theory. Since we want to apply continuum 
perturbation theory, we need additionally 

A</i. (1.4) 

A different renormalization method exploits finite-size effects and it is based on the 
observation due to Symanzik ^ that renormalizability is not spoiled when a field theory 
is considered on a finite space-time manifold. We can define a renormalization scheme 
by using the lattice size L (instead of the external momentum employed in the preceding 
scheme) as the relevant energy scale. Finite-size scaling schemes have been much studied 



in the last years ||8Hl4|. They have been successfully applied to the computation of the 
A-parameter, of the light-quark masses and of the renormalized axial current in quenched 
QCD. With respect to the previous method, there is the advantage that the requirements 



( |1.3| ) reduce to the weaker ones 

L > a, Aa < 1, /la -C 1. (1.5) 

Recently, a new method for the determination of renormalized composite operators in 
asymptotically free theories has been proposed in |1T5|-[18|. It is based on the Operator 



Product Expansion (OPE), which has been postulated for the first time by Wilson [19 



thirty years ago. It has been proved in perturbation theory by Zimmermann and it 
is widely thought to hold beyond perturbation theory. This new procedure shares many 
features of the infinite-volume schemes and in fact it is defined in infinite volume, but 
it has the advantage of being more direct. It provides the matrix elements directly in 
continuum schemes and therefore, it could be the method of choice in those cases in which 
the renormalized operators are obtained as combinations of many different lattice operators. 
Moreover, it avoids the evaluation of products of local operators at coincident points. This 
fact reduces lattice artifacts and allows a simpler implementation of operator improvement. 

We proceed now to describe the general context to which this method applies. Let us 
consider, for instance, the following simple example of OPE: 

A{x)B{-x) Co{x)O{0) + . . . (1.6) 

where the dots . . . indicate terms of higher order in x^, corresponding to operators of higher 
canonical dimension. Let us suppose that we know how to define the renormalized operators 
A and B nonperturbatively. For instance, this is the case if they are conserved currents: If 
the lattice does not break the corresponding symmetry, one can choose a discretization in 
which they are exactly conserved and therefore, do not renormalize. 

Let us suppose that we are interested in computing some hadronic matrix element of 
O which we shall denote by {hi\0\h2). The proposed procedure works as follows: 

1. The Wilson coefficient is calculated using RG-improved perturbation theory. Any 
renormalization scheme can be used in this step, but we will usually consider the MS 
scheme. Moreover, to avoid spurious scale dependences, we consider RG-invariant 
(RGI) Wilson coefficients and correspondingly RGI operators. We indicate by Cq (x) 
the resulting /-loop approximation for the Wilson coefficient. 
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2. The matrix element {hi\A{x)B{—x)\h2) is computed in a numerical simulation for a 
properly chosen range of x, say for < x^ < R^. This step gives a function G'^e(x). 

3. Finally, Gji^Bi^) is fitted using the form Cq{x) ■ O and keeping O as the parameter 
of the fit. 

The result O will be a function of Aa, L/a, and of the considered range [r^,i?^]. The 
physically interesting quantity^ is then 

0";?^,)= lim lim lim 3. (1.7) 

_RA,rA^O Aa,a/r,a/_R-+0 L/a-*oo 

Note that the limits cannot be interchanged. First, we should consider the infinite- volume 
limit, then the scaling (continuum) limit at fixed -RA, and rA, and finally take the short- 
distance limit. The outcome of this step,^ 0™"^*, is identified with the matrix element we 
are looking for, {hi\0\h2), renormalized in the same scheme in which we computed the 
Wilson coefficient Cq{x). 

In practice, one works at finite values of a and L. However, one may still hope that 
there exists a window of values of such that a ^ |a;| ^ A^^ ^ L, so that |x| is small 
enough to allow the application of perturbative OPE and still large enough to make lattice 
artifacts small. The volume must of course be large enough to avoid size effects. Notice 
that the condition a -C |x| <^ A~^ replaces here the condition a <C -C A~^, fi being 
the momentum at which renormalization is performed, present in standard infinite- volume 
schemes. In this respect the OPE method is not different from the standard renormalization 
methods. 



The OPE has already been used in lattice simulations. For instance, in Refs. p^-pG 
the authors consider quark and gluon two-point functions in QCD. Their aim is to compute 
the (chiral and gluon) condensates. They use an "OPE-inspired" phenomenological fitting 
form, considering the condensates as fitting parameters. In all these and similar studies 
vacuum expectation values of short- distance products are considered. The possibility of 
varying the external states, investigated in this paper, provides a stronger check of the 
OPE and could be a fruitful resource. 

An approach which allows to bypass the perturbative computation of the Wilson coef- 
ficient (with the related limitations) has been proposed in Refs. p7|,P^ . The idea has been 



tentatively applied to the evaluation of the momenta of the nucleon structure functions. 
This method allows to study the structure functions at energy scales <C 1/a^. 

In this paper we perform a feasibility study of the procedure outlined above in the two- 
dimensional nonlinear cr-model. This model is quite interesting since it shares with QCD 
the property of being asymptotically free with a nonperturbatively generated mass gap 
||29|^I1 . It has also been extensively studied in perturbation theory, both in the continuum 



|3^-p4[ and on the lattice p5|-^. Extensive simulations have been done to check whether 



^One could also consider the OPE in a finite physical volume. In this case, one would take the infinite- 
volume limit and the continuum limit together, keeping the physical volume fixed. 

^ Notice the close resemblance with the so-called "point-splitting" renormalization procedure proposed 



a long time ago by Dirac |21 
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the theory shows the correct behaviour predicted by the perturbative RG |38|,^3|- ^j and 
to test the exact predictions for the mass gap obtained using the thermodynamic Bethe 
ansatz Another advantage of this model is the availabihty of a very efficient 

algorithm ^,0,^, which shows no critical slowing down (for a general discussion see 
Ref. |5^). Therefore, we can work on very large lattices keeping the size effects completely 
under control. 

We consider the OPE of products of Noether currents, a particularly simple case since 
the current is exactly conserved and therefore does not need any renormalization. We 
discuss in detail all possible sources of systematic errors, by considering several different 
cases. 

The paper is organized as follows. In Sec. ^we define the model and the basic notations, 
and in Sec. ^ we give the OPE for the product of the Noether currents that will be studied 
numerically. Sec. | presents a detailed discussion of our numerical results. Our conclusions 
are presented in Sec. ^. In App. |^ and ^ we give our perturbative results, while in App. 
we critically discuss different methods for the determination of the running coupling 
constant. 

Preliminary results of this work have been presented at the Lattice conferences in 1998 
and 1999 



2 The lattice model 

In this paper we consider the nearest-neighbor lattice nonlinear cr-model in two dimensions. 
The fields are unit-length spins rr^ G -R^ and the action is given by 

where (9^/^ = {fx-^a^ — /a:)/o, and is the bare lattice coupling constant (often one 
introduces /? = l/gLj- The partition function is given by 



Z = j exp(-S^-") \[d^(T^5{al - 1) . (2.2) 



As usual, we introduce the lattice A-parameter 



1 



h}^'\a,gL) = - (/3o^7L)-^^/^«exp 
a 



PoQL Jo KP'^Ht) Pot' P'otJ\ 
1 /(iV-2)^?A-^/(^-^) f 2n 



a 



where P^^^^igi) is the lattice /5-function which can be expanded as 



k=0 



5 



The coefficient f3]^^^ are known up to four loops, i.e. for k < 3, see ^5, 37, 41]. The ffist 
two coefficients are universal, and therefore, we have not added the superscript. 
The mass gap of the model is related to the A-parameter by 



(2.5) 



The constant Cat is a nonperturbative quantity, and for generic models, it is usually un- 
known. For the nonlinear cr-model, an explicit expression has been obtained using the 
thermodynamic Bethe ansatz HO,BU|. Explicitly 



N 



l/{N-2) 



r(i + ]vb) 



25/2 



exp 



TT 



2{N -2) 



(2.6) 



In this paper we will consider the OPE of products of the Noether 0{N) current, which is 
defined by 

^ {<d,al - . (2.7) 



■L,ab 



9l 



Because of the 0{N) invariance of the theory, this current satisfies the Ward identity 



—a' - —a- 



(2.8) 

holds exactly on 



where C is a generic operator and {d~f)x = {fx — fx~afj.)/ci- Since 
the lattice, j^'"'' does not need to be renormalized. 

A second operator we will be interested in is the energy- momentum tensor. In the 
continuum it is given by 



T, 



1 
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1 2' 
df,(T ■ dycr - -6^u {da) 



(2.9) 



The energy-momentum tensor is the Noether current of the translations, and therefore it is 
exactly conserved. On the lattice translation invariance is lost, and thus, there is no lattice 
operator whose naive continuum limit corresponds to (|2.9|) and that is exactly conserved. 
However, as shown in Ref. [^], it is possible to define an operator which is conserved in 
the continuum limit.0 Explicitly one defines 



platt 

HUjX 



yL,i2,0)rpL I 7^^,(2,0) 
^TT ^ tJ.u,x "T ^T2 



^ 5,^{d,a)l + Z^-t'^-5,^{dar, 



9l 



9l 



(2.10) 



where {d^f)x = {fx+an — /x-a^t)/(2a), T^,^,^ is the naive lattice energy-momentum tensor. 



rpLi 



{d^(T)x ■ {dy(T)x - -5^,y{d(T)l 



(2.11) 



The general ideas have been presented in [|j , fohowing closely the method used in for the definition 
of the axial current. 
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and Z< 



L,{2,0) ryL,{2fi) 



Z^'^^'""' are renormalization constants. In Eq. ( 2.101) we have not 



;l,(2,o) 

considered terms that are not 0{N) invariant and that vanish because of the equations 
of motion.^ We have also discarded the contribution proportional to the identity operator 
(diverging as 1/a^), since it does not contribute to connected correlation functions. The 
renormalization constants have been computed to one-loop order in We have found a 
small error in one of the expressions. The correct results are reported in App. |A.2 . 



3 OPE for the 0{N) Noether currents 

In this Section we report the general form of the OPE for the products of the 0{N) currents. 
The coefficients are computed in perturbation theory, the explicit expressions to one-loop 
order being reported in App. ^ As is well known, a perturbative expansion cannot be 
defined directly for the model ( p.l| ) because of infrared singularities. A standard way out 
consists in adding a magnetic field that breaks the 0{N) invariance of the theory and plays 
the role of infrared regulator. One considers therefore 



S 



latt 



ha' 



291 



9l 



(3.1) 



A perturbative expansion is obtained setting cr^. = (tt^;, a/1 — tt^), and expanding in powers 
of 7r„. When the integration measure is kept into account, the total action reads: 



N 



(3.2) 



The calculation of the OPE is fairly standard and we report below the results in the 
continuum MS scheme and on the lattice for the cases that will be studied numerically. 



3.1 The continuum MS scheme. 

We begin by considering the OPE of the scalar product of two currents. The general form 
of the OPE is dictated by 0{N) and rotational invariance and by power counting. The 
result can be written as follows: 



a,b 



WAx) + 



rf> rf rf rf> ry rf> A I rf rf A 

^ ^ W2(x) + ^ ^W3(x) + 



[X 



2^2 



x^ 



^Similar terms appear also in QCD. Indeed, in this case, one should also consider, beside gauge-invariant 
operators, operators which are BRS variations and operators that are proportional to the equations of 
motion, see 0]. 
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+ 
+ 



Wrix) + 6ppWsix 



x^ 



^Hm5(0) + 



+ ^Wo,f,p{x)-l, 
x^ g 



(3.3) 



where Wi{x) and Wo^fj^p{x) are functions of x, of the MS couphng g, and of the renormal- 
ization scale fi. Exphcit one-loop expressions are reported in App. |B.1| . With [ ■ we 
indicate the MS renormalized operator, while T,ycr{x) is the energy-momentum tensor 
Notice the appearance of the operator a{x), 



aix] 



X 



(3.4) 



where the subscript B indicates bare quantities. This operator is not 0{N) invariant and 
its presence is due to the breaking of the 0{N) invariance by the magnetic term. However, 
using the equations of motion, a{x) can be rewritten as 



6S_ 
6n 



(3.5) 



Therefore, on-shell and in the limit /i 0, we recover an 0(A^)-invariant expansion. Ex- 
phcitly: 



2^^^{X)■Jp{ 



-XI 



^ apX yX cr 



Wi{x) + 



p-^p<^l''^C 



[X 



2^2 



W^ix) + 



W.(x)+ 



WAx) 



+ 



X 

'^Wl,{x) + 5,pWl,{x) 



x^ 



^l[(9.)^]-^(0) + 



with 



+\wo,pp{x)-l, 
x^ g 



W^(x) = W.(x) - Wvix) 



W'(x) = WAx) - WAx) 



The Wilson coefficients satisfy the following RG equations: 

Wo,pp{x;g,iJ,) = 0, 

Wi{x;g,iJ,) = 0, i = l,.. 

W-{x;g,fi) = 0, i = 5,6. 



^T.^f^^'^d-g-^ 
9 .9 (3{g) 

'T,^^^'^d-g-^ 



9/i dg g dg\g 



(3.6) 
(3.7) 

(3.8) 
(3.9) 
(3.10) 



where I3{g) is the MS /^-function. 
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We also consider the OPE of the antisymmetric product of the currents. Neglecting 
terms of order 0{x\og^ x), we have^ 



[X 



2\2 



X^ 



X" 



H^,J.p - 5,MUr{x)- (0) - ^^^^(O)] 

The coefficients Ui{x) and Uoi{x) satisfy the RG equations: 

+ ^ 



(3.11) 



(3.12) 



3.2 Lattice 



On the lattice one can write expansions completely analogous to those holding in the con- 
tinuum. The only difference is the appearance of terms that are cubic but not rotationally 
invariant. For the scalar product of two currents we obtain 



-7 • 7 



^ apX yX o 



'^^■^-f^W,\x) + ^I^^f^W^ix) + tlEtJ^^i£l^_:i£t^w^{x) + 



X uXjyS pfj ~\~ X pXfjSfj 



x^ 



[X^ 



X^ 



+ 



-\- 



9l 



.2 + (r,2\2 



ryt ry rf rf /y ryt A I rp rp A 



X^ 



[X^ 



X^ 



^ pumper + ^P°"^P'" ^^^^^ 



91 



+ 

+ 



''''''''W,'^{x) + 5,,W,\x)]^idcT)l + 



x" 

'^w,Hx) + 5,M{x) 



x^ 



''pp''" 8 



9l 
9l 



+ \w,^,,,{x)-l. 
X 9l 



(3.13) 



Here Tj^.^, is the lattice energy- momentum tensor defined in Eq. ( p.ll| ), and is the 
lattice analogue of the quantity defined in Eq. 



a. 



N 



h + d'a 



9l 



(3.14) 



Note that one could also add a contribution proportional to [x^x^S, 



up 



x^Xa^ ^(3) / dcxj^ (0) ^f3ja (^) ' Howevcr, in two dimensions, {x^^Xa^up 

(a ^ (3)) 

x'^{Sfj_aSuf3 — 5va5^p), and thus this term is equivalent to that proportional to Ui(x). 
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where d'^ = d~dfj,. Note the last term in Eq. ( 3.14|) that is not present in the continuum 
and which is due to the measure term in the functional integral. Using the lattice equations 
of motion, one can show that 



f Clatt „ 



(3.15) 



Therefore, on-shell and in the limit h 0, it is possible to get rid of the terms proportional 
to as we did before. 

In the antisymmetric sector we can write 



E\-L,ac-L,bc 
^Jfj.,x Ju,-x 



■L,6c -L.ac 
J/i,a; Jv,—x 



rf rtft np 



9l 



■L,ab 



(3.16) 



3.3 RG-improved Wilson coefficients 



We want now to resum the Wilson coefficients using the RG. We assume that the Wilson 
coefficient W{^x; g) satisfies the equation 



W{^ix-g) = 0, 



(3.17) 



where P{g) is the (scheme- dependent) /3-function and 7^(5') is related to the anomalous 
dimensions of the operators appearing in the OPE. Such an equation is not valid in general, 
since for generic operators 'y^{g) is a matrix. However, in the cases we are interested in. 



cf. Sec. |3]^ and |3^ , we can consider 711/(5') as a scalar function of g and thus apply Eq. 
The general solution of Eq. ( ^.17| ) is easily found. If we define a running coupling g{g,t) 

by 



exp 



(3{x) 



then 



^^7(^r 



W{iJ,x; g) = W{1; g{g, ^ix)) exp 

It is useful to rewrite this equation as 

W{iix]g) = U{g)WRGi{g{g,fxx)) 



(3.18) 



(3.19) 



(3.20) 
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where 



WRGI{9{9,^^x)) = W{l■g{g,^lx)) [g{g , ^lx)]^^ ' exp 



dx 



Uig) 



— dx 



+ 



w 
7o 



(3.21) 



where we have expanded 7^(5') = 7^(7 + 0{g'^ 
definition of the operator, introducing 



The factor U (g) may be included in the 



Orgi= g''"/^'exp 



— dx 



(3{x) l3o 



(3.22) 



where 7^(5^) = 7^(7 + 0{g^) is the anomalous dimension of O. The new operator Orgi is 
called RG invariant, in the sense that it satisfies a RG equation of the form 



On 



On, 



0, 



(3.23) 



where ^'o\gi n-point irreducible correlation function with one insertion of O. In 

other words, with this new definition, we have eliminated the running due to the anomalous 
dimension of the operator. Note, however, that, in spite of the name, this definition is not 
scheme independent. 

Eq. ( p. 181 ) can be written as 



K{l,g{g,^ix)) = Ain,g)x, 



(3.24) 



which shows that g is a function of A(yU, g)x only, i.e. g{g, fix) = g{Ax). Moreover, defining 
z = — log Ax, we can expand 'g{Ax) in inverse powers of z. Explicitly, we obtain 



(3.25) 



Since the (3 function is known to four-loop order, both on the lattice and in the MS scheme, 
two additional terms can be trivially added to this expansion, neglecting terms of order 
\og^ z/z^. 

Finally, we wish to point out a further interpretation of g. Consider for instance the MS 
scheme and correspondingly the MS coupling gjfg which is a function of m//i, where m is 
the mass gap. Then, it is easy to show using the RG equations that g{g, fix) = gj^ijnx). In 
other words, g{g,^x) is nothing but the coupling constant computed at the scale /i = 1/x. 

Finally, notice that the results presented in this Section can be applied to any renor- 
malization scheme. In the following, we will use them to resum continuum MS predictions 
as well as lattice perturbative expressions. 
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4 Numerical results 



In this Section we present our numerical computations meant to verify the feasibihty of the 
nonperturbative renormahzation method presented in the Introduction. 

We consider matrix elements of the product of two Noether currents between one- 
particle states. This case is particularly simple, because the currents are exactly conserved 
and the matrix elements of the operators that appear in the OPE expansion can be com- 
puted exactly in most of the cases. 

The typical procedure we adopt is the following: 

1. We compute a matrix element {l\A{x)B{y)\2) measuring a suitable lattice correlation 
function. 

2. We compute the matrix element (1|(9|2) appearing in the OPE expansion either 
exactly (as is possible in most of the cases) or numerically by means of some different 
numerical nonperturbative technique. 

3. We divide {l\Aix)B{y)\2) by the OPE prediction J^o^oix - y){l\0\2). 

The goal is to see if there is a window of values of \x — y\ in which the OPE works, i.e. the 
result of step 3 is 1 independently of |a: — For the cases we will consider here, the OPE 
gives an accurate description (at the level of 5-10%) of correlation functions for distances 
2^^ — 1/1^^ = (from now on we set a = 1). This result is quite encouraging for 
future applications of this method. 



4.1 The observables 

We have simulated the 0(3) a-model with action (|2.1| ) using a Swendsen-Wang cluster 



algorithm with Wolff embedding ||43| , pl|-|53| . We did not try to optimize the updating 
procedure: Most of the CPU time was employed in evaluating the relevant observables 
(four-point functions) on the spin configurations of the ensemble. In order to estimate the 
scaling corrections, we simulated three different lattices, of size T x L, using in all cases 
periodic boundary conditions: 

(A) . Lattice of size 128 x 64 with = 1/1.40. 

(B) . Lattice of size 256 x 128 with gi = 1/1.54. 

(C) . Lattice of size 512 x 256 with gi = 1/1.66. 

The algorithm is extremely efficient — the dynamic critical exponent z is approximately — 
and the autocorrelation time is very small. We performed a preliminary study in order to 
determine how many iterations are needed to obtain independent configurations. For this 
purpose we measured the normalized autocorrelation function 



^-3 Ef=iiOi-oy 



m = ."--'^m':"^' ' (4-1) 
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for different observables O. Here A/" is the number of Monte Carlo iterations, Oi is the 
value of O at the i-th iteration, and O the sample mean of O. In Fig. |I] we report A{t) 
for the observable 

Od = -^j^'^^io'x ■ f^x+v + (^x ■ (^x+w), (4.2) 

x 

where v = (d, 0) and w = (0, d). We have considered d = 1 — a short-distance observable — 
and d ~ ^''^^ = m~^. As expected, local observables have a slower dynamics — this is due 
to the fact that the dynamics is nonlocal — than long-distance ones. In any case, for all 
observables the autocorrelation function shows a fast decay: indeed, for c? = 1 we have 
^exp ^ j^g^ while for d ~ ^'^^p we have t'^^p ^ 5. Moreover, A{j) is independent of the 
lattice used, confirming the fact that z ~ (as we shall report below, all lattices have the 
same L/^'^^p). Since the measurement of the observables is quite CPU-time consuming, 
we computed the relevant correlation functions only every 15 iterations, which should be 
enough to make the measurements independent. Nonetheless, most of the CPU time is 
employed in evaluating the observables on each spin configuration of the ensemble. We 
computed the product j{0)j{x) of two Noether currents for distances x smaller than some 
fixed fraction of the correlation length: |x| < k^. If the physical size L/C, of the lattice 
is kept constant (as we did) we expect the CPU time to scale as L^. The CPU time per 
iteration turns out to be roughly independent of the particular product considered. As an 
example we give the CPU time per measurement for the simulation in which we compute the 
antisymmetric product of two currents between states with opposite momentum. For the 
three different lattices, on an SGI Origin2000, we have: ri28x64 ~ 5.4 sec, r256xi28 ~ 71 sec, 

T512X256 ~ 1100 sec. 

We measure several different observables. First, we measure the two-point function 
C{p]t) (here and in the following the "temporal" direction is the first one, of extent T) 

1 ^ ^ _ 

C{p;t)^—Yl E e^^(---)(^T,,.,-^,+,,.,). (4.3) 

tl=l Xl,X2 = l 



We computed the correlation function C{p;t) on the lattices |(B)| and |(C)| for momenta 



p = 2TTn/L, n = ... 3 and times separations < t < 100; on lattice |(A)| we considered 
the same set of momenta and time separations < t < 40. The number of independent 
configurations we generated is: Nconf — 6 ■ 10^ for lattice |(A)| ; Nconi = 590000 for lattice 



fBl ; A'conf = 180000 for lattice [C] and p = 0; finally A^conf = 139000 for lattice [C] and 
p^O. 

A check of our simulation is provided by the results of Ref. [^, who computed, among 



other things, the mass gap for lattices |(A)| and |(B)| . For the exponential correlation length 
^exp _ ^-1 obtain 

f = 6.878(2), 13.638(10), 27.054(25), (4.4) 



for lattices |(A)| , |( B )| , |(C)| respectively. They are in good agreement with the results of 



Ref. 0: they obtain ^^^p = 6.883(3) and ^^^p = 13.632(6) for the first two lattices. The 



three lattices we simulate have approximately the same physical size, mL ~ 9, which is 
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p 


O = (9oO-)2 


O = [doo- ■ ditr) 


O = {d^crf 




2n/L 
Atx/L 
Qti/L 


34.619(25) 
34.707(18) 
34.735(18) 
34.741(26) 


0.00053(48) 
0.03065(26) 
0.06080(39) 
0.09026(64) 


34.663(25) 
34.776(18) 
34.857(18) 
34.923(26) 



Table 1: Estimates of Y2a^cfiP^P'^ different operators measured on lattice |(B)| . For 

{docr ■ dicr) we report the imaginary part, the real part being zero. The matrix element of 
the other two operators is real. 



p 


O = {d,af 


O = {doo- ■ di(T) 


O = (9iO-)2 


27r/L 
A-n/L 
Qti/L 


0.25434(60) 
0.25153(85) 
0.25597(128) 


0.01452(47) 
0.02824(61) 
0.04024(96) 


0.20241(68) 
0.18394(94) 
0.16643(140) 



Table 2: Estimates of Yla^o^iP^ ^0) for different operators measured on lattice |( B )| . We 
report here the real part for (Sqct)^ and (9i<t)^, and the imaginary part for {docr ■ dicr). 



large enough to make finite-size effects much smaller than our statistical errors. This is 



confirmed by the analysis of Ref. [^7 



In order to verify the OPE, we need the values of the matrix elements which appear 
in the r.h.s. of Eq. (|1.6|) . Matrix elements of lattice operators can be computed from 
properly defined three-point correlation functions. If Ot^x is a lattice operator, we define 
the correlation function 

T L L 

to = l Xo = l Xl,X2 = l 

and the corresponding normalized correlation 

The function C^{p,q;2t) has a finite limit for t — > oo. In this work we will only need 
the matrix elements of the naive lattice energy-momentum tensor ( |2.11| ). For this reason. 



we have computed C^{p,p]2t) with O = df^cr ■ dpcr. Such a correlation function has 
been computed on lattice |(B)| , using Nconi = 320000 configurations, for t = 5, . . . , 20 and 
p = 2mi/L, n = 0, ...,3. For these observables CQ{p,p;2t) is independent of t, within 
the statistical errors, already at t = 5. The results obtained for t = 5 are reported 
in Table For O = Oqct ■ dicr, statistical errors are dominated by the error on the 
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evaluation of (^^(p, p;2t). On the other hand, for O = {df^crY, the statistical error of 
the numerator in Eq. ([4.6| ) is roughly equal to that of the denominator. The reason is 
clear: since in the continuum limit ((9^cr)^ is proportional to the identity operator, we are 
computing essentially (up to O(a^) corrections) the same quantity in the numerator and in 
the denominator with approximately the same statistics. The reported errors on the ratios 
are obtained using the independent error formula. For ((9^cr)^ smaller error bars could have 
been obtained by taking into account the statistical correlations between numerator and 
denominator. 

We also measured C^{0, p; 2t) for the same operators on lattice |(B)| , using Nconi = 62000 
independent configurations, for p = 2mT/L with n = 1, . . . , 3 and t = 5, . . . , 10. The 
normalized three-point function shows a plateau for t >10 when /z = z/ and for t > 6 when 
/i 7^ z/. The results obtained for t = 10 are reported in Tab. |^. In this case the statistical 
errors are dominated by the uncertainty on C^{0,p; 2t). 

In this paper we study the OPE of the product of two currents. Therefore, we have 
computed the following correlation functions: 

G^^\t,x;p,q;2Q ^ ^ ^ <(^t(o,o) " ■ (4-7) 

Xl,X2 

G^;Xt,x;p,q;2Q ^ Y.H:m&.) ' 

xi,X2 abc 

where j^'^* is the lattice Noether current defined in Eq. ( ^^7]) . Of course, we averaged over 
lattice translations. 

In particular we have measured: 

a) G'^^\t,x;p,p]2ts) using Nconi — 1-3 ■ 10^ configurations on lattice |(A)| and using 
Nconi = 58350 independent configurations on lattice |(B)| . 



b) ImG'('*)(t, X] p, 0; 2ts) using Nconf = 112000 independent configurations on lattice [[K 



c) ReG^ii{t,x;p, —p;2ts) using Ai'conf — 2.4 ■ 10^ configurations on lattice |(A)| , A^conf = 
69500 independent configurations on lattice |(B)| , and iVconf = 31550 configurations on 
lattice [U]]. 



d) ImG[,^(t , x; p, 0; 2ts) using A^conf = 41750 independent configurations on lattice |(B) 



In all cases we consider p = 27m/ L, n = 1,2, 3; = 7, 8, 9 and |t| < 5, |x| < 5 on lattice |(A) 



ts = 10, 11, 12 and \t\ < 8, |x| < 8 on lattice [BJ, and = 20, 23, 26 and |t| < 12, |x| < 12 
on lattice |((J) . 



Using the four-point correlation function determined above, we constructed the normal- 
ized ratios 

^(.),, _ _ „ , X _ G^-\t,x-p,q;2t,) 

G^ '[t,x;p,q;2ts) = — , (4.9) 

y/C{p;2QC{q;2tS 

which have a finite limit for tg — > oo. We verified that G^'^ {t, x; p, q; 2ts) is independent of tg 
in the range considered, and thus we have taken the estimate corresponding to the lowest 
considered value of ts as an estimate of G^'\t, x;p,q; oo). 
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In the paper we will usually consider averages over two-dimensional rotations, i.e., given 
a function f{zt,Zx), we consider 

7(^) ^ %^i/if)^, Q^^,)^e(\z\-r + l]9(r + l-\z\], (4.10) 



with \z\ = >y zf + z'^ and r = n + 1/2, n integer. 
4.2 One-particle states 

In the conventional picture the lowest states of the model are one-particle states transform- 
ing as 0{N) vectors. On a lattice of finite spatial extent L, we normalize the states and 
the fields as follows: 

{p,a\q,b) = 2a;(p)L<5%,5, (4.11) 



{P,a\a'M = s^^5^'e^-^\ (4.12) 

The function uj{j)), which is the energy of the state |p, a), and the field renormalization 
Z(j)) can be determined from the large- |t| behavior of the two-point function C{p;t): 

C(p; t) ~ ^^e-"®* for t > 1. (4.13) 

In the continuum (scaling) limit we have Z(j)) = Z independent of p and uj(j)) = \/p^ + m?. 
In Table ^ we report our results for the three lattices. In order to evaluate a;(p) and Z{p) 
we determined effective values at time t by solving the equations 

C{p;t + l) _ cosh[^eff(P,t)(t+l-T/2)] 

C{p-t) C0Sh[^efT(P,t)(t-T/2)] ' ^ • ^ 

2uJcs{p,t) 

Then we looked for a plateau in the plot of co'efr(P, t) and Zef[{p, t) versus t. Both functions 
become independent of t for t > C,- The values reported in Table ^ correspond to one 
particular value of t of order ^'^''P: t = 8, 16, 20 respectively for lattice |(A)|, |(Bj|, |((J)|. For 



the two largest lattices there is evidence of scaling at the error-bar level. Instead, for lattice 
(A)| there are tiny scaling violations, that are however so small (at most 1%) that we can 



neglect them in the following discussion. 

One can also investigate asymptotic scahng, i.e. the dependence of and Z on 
gi. The dependence of 00 (p) can be determined from Eq. ( p.5| ). There exists also an exact 
prediction for Z, including the nonperturbative constant [^,^. However, as is well known, 
lattice perturbation theory is not predictive at these values of the correlation length, and 
indeed, the perturbative four-loop predictions show large discrepancies compared to the 
numerical data. The agreement is instead quite good |^,|5^ if one uses the improved 
coupling gE defined in Eq. 
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lattice 


(A) 


lattice (B) 


lattice 


(C) 


p 


uj{p) 


Z{P) 


uj{p) 


m 


uj{p) 


Z{p) 




2tt/L 
Att/L 
67r/L 


0.145393(40) 
0.175380(40) 
0.243657(74) 
0.326327(192) 


1.6593(8) 
1.6582(8) 
1.6510(13) 
1.6378(32) 


0.073327(55) 
0.088244(67) 
0.12263(15) 
0.16415(32) 


1.3563(18) 
1.3557(21) 
1.3591(42) 
1.3499(86) 


0.036963(34) 
0.044348(35) 
0.061456(66) 
0.082494(109) 


1.1295(14) 
1.1284(13) 
1.1292(21) 
1.1336(35) 



Table 3: The one-particle spectrum and the field normalization for lattices |(A)| , |(B)| , |(C) 



4.3 Nonperturbative renormalization constant for the lattice energy- 
momentum tensor 

In this Section we want to compute nonperturbatively the renormalization constant of the 
lattice energy- momentum tensor ( p.lO| ). In general, given an operator O on the lattice we 
define its matrix element by 



{p.a\d%h) = N^A^{p)u{q) lim C^\p,q;2t). 



(4.16) 



For d^cr ■ d^cr the matrix elements can be obtained from the results given in Tables |l| and ^ 
The matrix elements of {d^crY are dominated by the mixing with the identity operator and 
thus, in order to define the renormalized operator for fi = u, we should perform a nonper- 
turbative subtraction of the large l/a^ term, which is practically impossible.^ Therefore, 
we only compute the renormalized operator for fi u, which amounts to determining the 



constant Z, 



L,(2,0) 
TT 



This constant is obtained by requiring 

(p,a|T^f*|p,6) 



(4.17) 



In practice, we first compute an effective (momentum-dependent) renormalization constant 

2puj{p) ql 



cm 



Im (p, a\{docr ■ dicr) \p, a) 



(4.18) 



which, in the continuum limit, becomes independent of p and converges to 



Using 

the data of Table |, on lattice [B], we obtain C{p) = 1-040(8), 1.048(6), 1.059(7) for 
p = 27t/L, 47r/L, Gir/L respectively. Clearly the scaling corrections are small, and thus we 



can estimate Z, 



L(2,0) 
TT 



1.05(2) on this lattice. This compares very well with the result of 



one-loop lattice perturbation theory given in Eq. ( |0^ ), which yields Z^J?'"^ ~ 1.044357, 



and with the result of "improved" (sometimes called "boosted") perturbation theory in 



terms of the coupling qe, cf. Eq. ( |C7| ), Z, 
on qe). 



L(2,0) 
TT 



1.052471(3) (the error is due to the error 



^In general we expect {p, a\{docr)^\p,b) — 2AL\/j? 



m? + B - Cii^ + m2) and {p,a\{di(Tf\p,b) = 
2 + B + Cp^ . The quantity we are interested in is C. However, from Table ^ we immediately 
realize that much smaller errors are required to really observe the momentum dependence of the matrix 



2AL^/^ + TO 



elements and thus to determine the constant C . 
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4.4 OPE for the scalar product of currents 

In this Section we consider the product Jq oil averaged over rotations. Using Eq. 
we have in the continuum scheme 



1 



^ab ^ab 
JO,OJl,x 



'-W2{r/2) + W3{r/2) + W^{rl2) 



[Toi(O)] 



MS 



(4.19) 



All other contributions vanish after the angular average. Using the results of App. |B.l 
have at one loop in the MS scheme 



we 



■ab 'ab 

:jo,oji,x 



1 - 



N-2 
~27 



9 (log ( 



/ir 
Y 



+ 7- 



4 



9 



[Toi(O)] 



MS ' 



(4.20) 



where /z is the renormalization scale and 7 Euler's constant. We will not use this form of 
the OPE expansion, but instead the RG-improved Wilson coefficients. Thus, cf. Sec. 
we write 

(4.21) 
(4.22) 



■jSfljt = WnoM^Msr)) [^01(0)]^^, 



where 



5(iV-2) 



in 



and g{Ajj^) is the running coupling constant defined by| 



Using the perturbative expression (|3.25|) , we can also rewrite (|4.22|) as 

N-2 



WRGiir) = (3oZ + 



Pi , ^ 5(iV-2) 
— logz + 



2tt 



-z + 



1 , 5 AT -2 

— log ;z + ^ - 

2n ^ Stt 



(4.23) 



(4.24) 



where z = — log(A^jgre'''/2). 

We will also use the lattice Wilson coefficients defined in Sec. |3^ . Using the one-loop 
results of App. p.3| and the general expressions reported in Sec. |3]^, proceeding as before, 
we obtain the prediction 



Jofjtf = m9L)Wliaii9Li^'^''r))T,\. 



(4.25) 



where, at one loop. 



W^Gli9L) = =- 

9l 

U{9l) = 1 + 



1 + 



5A^-2 



in 



1]9l 



1 1 

n 4 



(4.26) 
(4.27) 



^This equation should be intended as follows: Aj^(l, g) = Aj^(/i, g)/ 11 in the right-hand side is assumed 
given, while for Aj^(l,5) we intend its expression in terms of the /3-function, Eq. 
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and (yf^(A^'^**r) is the running coupling constant defined by 

A'^"(1,^J= A'^"(l,^7i)re^v^, (4.28) 

where A^^** is the lattice A-parameter ( p.3|) . 

Finally we shall test perturbation theory in the "improved" expansion parameter 
defined in Eq. (|C.7|) . The OPE becomes: 

1 



2 Joh Ji:: = U{9E)W^aii9E{^Er))T,l,, (4.29) 



where 

1 r /^(M - ')\ ^\ 1 

(4.30) 



WiciigE) = — 
9e 



^ , 5(iV-2) 1, 



[/(■) is the same as in Eq. ( [4.27| ), and g^i-^E^) is the running coupling constant defined by 

Ae{1,9e)= AE{l,9E)re^V8, (4.31) 

where A^; is the A-parameter (|C.9|) . 

We have tested the validity of the OPE by considering matrix elements between one- 
particle states. The matrix elements of the product of the currents]^ can be determined in 
terms of G^^\t, x;p,q; 2ts), since 

^ (P.c|jo_(o,o) ■ Ji,(t,x)l^,c) = N^Auj{p)u{q) lixn^ G^'\t, x; p, g; 2t,). (4.32) 



In Fig. we report]^ the angular average of Im G^'^\t, x; p, p; 2ts) for lattices |(A)| and |(B) 



here ts = 6, 10 for the two lattices respectively. In Figs. |, and ^ we compare these 
numerical data with the predictions of perturbation theory. 

In Fig. ^ we use continuum RG-improved perturbation theory in the MS scheme. In 
this case, the matrix element of the energy-momentum tensor is immediately computed: 
{p,a\TQi\p,b) = 2ipuj{p)6"'^. Therefore, we consider the ratio 



l Im(p,c|jo,(o,o)-ji,(,,,)b,c) 
""^"^ 2 2MmRGimMsr)) ' ^ ' 

which should approach 1 in the short-distance limit. In Fig. ^ we present several determi- 
nations of R{r) that differ in the way in which the running coupling constant ^(A^^r) and 
the MS coupling g are determined. 

There are several different methods that can be used to compute the MS coupling g and 
the strictly related ^(Aj^x). They are compared in detail in App. 0. It turns out that the 

^ Since the currents are exactly conserved, both on the lattice and in the continuum, there is no need 
to make a distinction between lattice and MS -renormalized operators. 



^° On lattice (B) ^^^t, x;p,p;2ts) has only been measured in I? = {(i, x) : \t\,\x\ < 8}. The points with 
r > 8 appearing in the figure correspond to "partial" angular averages, i.e. they have been obtained using 
Eq. (4.10) and restricting z to V. The same comment applies also to the subsequent figures. 
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finite-size scaling method proposed by Liischer 0,61] and what we call "the RG- improved 
perturbative method" are essentially equivalent, see, e.g.. Table § in App. |C[ Therefore, we 
can use either of them,[^ obtaining completely equivalent results. In Fig. ^ we have used 
the finite-size scaling method to be consistent with what we would do in QCD. Elsewhere, 
we have used the RG improved perturbative method because of its simplicity. 

The first step in the computation of R{r) is the determination of ^(Aj^r). This is 
obtained as follows: we fix fi/m, and using the finite-size scaling results reported in Table 



in App. y, we compute QMsifi). Then, using Eq. (|C.4|) we determine Aj^ at /-loops. 



Finally, ^'(Aj^jgr) is obtained either by solving Eq. ( |4.23|) , again using Eq. (|C.4|) for Aj;jg(l, g) 



appearing in the right-hand side, or by using Eq. ( |3.25|) . As we discussed in Sec. 

the final result should be independent of the chosen value of /i/m and therefore we can 

evaluate the systematic error on 'g{A-^x) by considering different values of /x/m. If we fix 



m//i = 0.00071, cf. first row of Table § in App. 0, on lattice |(B)| , we have gjA^^x) = 1.498 



2.206, 3.363 respectively for mx = 0.2, 0.5, 0.8, while if we fix m//i = 0.1033, cf. 10th 
row of Table § in App. 0, we have glAj^x) = 1.490, 2.185, 3.273 at the same distances. 
The dependence is tiny and, as expected, it increases for larger values of mx. In practice, 
it does not play any role, the main source of error being instead the truncation of the 
OPE coefficients. Notice that the independence on n is obvious if we use the RG-improved 
perturbative method. 

Let us now describe the various graphs appearing in Fig. In graphs (A) and (B) 
we fix /i = m/0. 00071 and then, using the results presented in Table ^, cf. first row, we 
obtain 5'Ms(/i) = 0.587016. We then compute Aj^(/i, (yfj^) using the four-loop expression 
with / = 4. Finally, the Wilson coefficient is given by ( [4.24 ): in graph (A) we use the 



leading term only, while in graph (B) we include the next-to-leading one. 

In graphs (C) and (D) we compute A^ as in (A) and (B), choosing a different scale, 
fi = m/0. 1033, cf. lO*'^ row of Table ^ Then we compute g{x) by solving numerically Eq. 
( [4.23| ) using the four-loop expression ( |C.4| ) with I = 4. Finally, the Wilson coefficient is 
obtained using Eq. ( [4.22| ). While in graph (C) we keep only the leading term l/g(x), in 



graph (D) we use the complete expression. 

Finally, graphs (E) and (F) are identical to graph (D) except that Aj^ and g{x) are 
computed using the two-loop expression ( |C.4|) with 1 = 2. The two graphs correspond to 
different choices of /x: /i = m/0. 00071 (graph (E)) and /i = m/0. 1033 (graph (F)). 

Comparing the different graphs, we immediately see that the choice of scale fi and the 
order of perturbation theory used for A^ (two loops or four loops) are not relevant with the 
present statistical errors. Much more important is the role of the Wilson coefficients. If one 
considers only the leading term (graphs (A) and (C)) there are indeed large discrepancies 
and in the present case one would obtain estimates of the matrix elements with an error 
of 50-100%. Inclusion of the next-to- leading term — this amounts to considering one- loop 
Wilson coefficients and two-loop anomalous dimensions — considerably improves the results, 
and now the discrepancy is of the order of the statistical errors, approximately 10%. For 

^^In QCD the RG-improved perturbative method cannot be used since no exact prediction for the mass 
gap exists. In this case the finite-size scahng method would be the method of choice. Note that for large 
values of the scale also the improved-coupling method [^,^ , which can be generalized to QCD |64j , works 
well, see App. ^. 
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the practical application of the method, it is important to have criteria for estimating the 
error on the results. It is evident that the flatness of the ratio of the matrix element by 
the OPE prediction is not, in this case, a good criterion: The points in graph (A) show a 
plateau — and for a quite large set of values of r — even if the result is wrong by a factor 
of two. However, this may just be a peculiarity of the case we consider, in which the r- 
dependence of the data and of the OPE coefficients is very weak. On the other hand, the 
comparison of the results obtained using different methods for determining g{x) seems to 
provide reasonable estimates of the error bars. For instance, if only the leading term of 
the Wilson coefficients were available, we could have obtained a reasonable estimate of the 
error by comparing graphs (A) and (C). 

In Fig. ^ we consider lattice RG-improved perturbation theory, computing 



^ ' 2U{g,)W^^,{-g,)lm{p,c\Tl^,\p,cy ^ ' > 

In graph (A) we use ql as an expansion parameter. We compute A^^^^{1, ql) using the 
value of the mass gap = 13.632(6) and Eq. ( |2.5| ), with the nonperturbative constant 
( p.6| ). Then, we determine 'gi^{A^^^^r) by solving numerically (|4.28|) and using for A^'***(l,^j^) 



appearing in the left-hand side its truncated four-loop expression (|2.3| ). Finally, we use Eq 



( [4.26| ) for the Wilson coefficient. The results are quite poor: there is a downward trend as 



a function of r and the data are far too low. Naive lattice perturbation theory is unable to 
provide a good description of the numerical data. 

The results improve significantly if we use the improved coupling qe'- The estimates 
obtained using this coupling, graphs (B), (C), (D), are not very different from those obtained 
using MS continuum perturbation theory. Graph (B) has been obtained exactly as graph 
(A), except that in this case we used qe as an expansion parameter. The A-parameter is 
defined in ( |G.9| ) , m/ A.e in ( |G.10| ) , and the relevant Wilson coefficient is given in Eq. ( |4.3CI|) . 



Graphs (C) and (D) are analogous to graph (B). The difference is in the determination of 
Ae{^i Qe)- We do not compute it nonperturbatively by using the mass gap but we determine 
it directly from the perturbative expression (|C.9D . In this case we use the perturbative 
expression truncated at four loops (C) and two loops (D). 

In Figs. 1^, ^ we have checked the validity of the OPE on lattice |( B )| . If one has in mind 



QCD applications this is quite a large lattice since ^''^^ ~ 14. For this reason, we have 
tried to understand if the nice agreement we have found survives on a smaller lattice, by 
repeating the computation on lattice |(A)| . The results are reported in Fig. ^. Graphs (A) 
and (B) should be compared with graph (D) of Fig. |^. In (A) we compute Aj^ from Eq. 
( ^75|) , with the nonperturbative constant ( p.6[) and = 6.878(3). Then, we compute 
'g{x) solving numerically Eq. ( [4.23|) using the four-loop expression ( |C.4| ) with I = 4. The 
Wilson coefficient is obtained from Eq. (|4.22 ). In (B) we repeat the same calculation as 



in Fig. ^ graph (D) at the scale /i = m/0. 00071. In (C) and (D) we repeat the calculation 
of graph (B) using the two-loop and the three-loop /^-function for the determination of the 
couphng g. In all graphs (B), (C), (D) we use the result = 6.878(3). 

Graphs (A) and (B) show a nice flat behaviour and for 2 < r < ^, the ratio R{r) 
is approximately 1 with 2-3% corrections (notice that the vertical scale in Figs. ^ and 
^ is different): The OPE works nicely even on this small lattice. (A) and (B) differ in 
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the method used in the determination of the MS couphng. As we explained in App. |y 
and it appears clearly from the two graphs, the effect is very small. Graph (C) — and to 
a lesser extent graph (D) — shows instead significant deviations: Clearly, ^(Aj^x) must be 
determined using four-loop perturbative expansions in order to reduce the systematic error 
to a few percent. Notice that such discrepancies are probably present also for lattice 



however, in this case, the statistical errors on R{r) are large — approximately 5-6% (for 
p = 2tt/L and p = Air/L) and 9% (for p = Gvr/L) — and thus they do not allow to observe 
this effect. 

As a further check we considered matrix elements between states of different momentum. 
In Fig. we report the angular average of Im G^'\t,x; 0,p; 20) for lattice [[B]]. In Fig. 
we compare the numerical data with the OPE prediction, by considering 



^ l Im(0,c|jo,(o,o)-j\^(,,,)b,c) 
) 2 iy«f?^(^)lm(0,c|r^f |p,c) ' ^ ^ ^ 



where T^t** is defined in Eq. ( p.lC ) and Z^'^^'^^ is computed in Sec. In Fig. |^ we report 



S{r) for lattice |(B)| . The Wilson coefficients are computed as in graph (A) of Fig. |^, using 
m^^ = 13.636(10). The numerical data are again well described by the OPE prediction for 
a quite large set of values of r. 

4.5 OPE for the antisymmetric product of currents 

In this Section we consider the antisymmetric product of two currents and compare our 
numerical results with the perturbative predictions. With respect to the previous case, we 
have here a better knowledge of the perturbative Wilson coefficients — some of them are 
known to two loops — and moreover we have an exact expression for the one-particle matrix 
elements of the current jf^^x- Indeed, we have [^|65[| : 



(P, c\j;,%Jq, d) = -z{p, + q,) G{k) (5-5^<^ - S'^'S"^) e'^^-") (4.36) 

where p ■ x = Pot + PiX, p^ = {i\/p^ + m'^,p), k = |a/(po - QoY + {Pi - QiY, and, for 
N = 3, 

"^'^ = 21^2-^' '«^) 

where the rapidity variable 6 is defined by A; = msinh0/2. 

We first consider the product {j^^o o)jt\t x) ~ i^i^{of))ii\tx)) for /i = z/ = 1 and x = 0. 
The OPE of such a product can be determined from Eq. ( p.ll|) . Using the results of App. 
B.2|, we have for t ^ 0, 



E 3um3um - itm^um = (^(Ams^)) 6o'(o) + ^^'lo)) , (4.38) 

where, at two loops, 

iV - 2 N -2 
WnGiig) = + -^9, (4.39) 
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and g{Ay^t) is defined in Eq. ( [4.23|) . We also consider the angular average of the product 
of the currents for fi = and u = 1. Using the results of App. |B.1| , we have for r — 



-T-aOf \ — -ac Abc Abc Aac ^ 

V J — Jo,{Ofi)Jl,{t,x) J0,{0,Q)Jl,{t,x) ^ 

' [dojt\0) + «'(0)] + + ] [%f (0) - «^(0)] . (4.40) 



167r ^ "-^^ ^ ' ^-^^ ^ 2g\ A-n 

Again, we have tested the validity of the OPE by considering matrix elements between 
one-particle states. The matrix elements of the product of the currents are obtained from 



E 

abc 



(P'«l4:t.4S,^-4o:i,.4:lJ^ = A^V^MM^ li^n G(t)(t,x;p,g;2t,). (4.41) 



In Fig. H we show a plot of Re C'fi (t, x; p, —p\ 2ts) obtained on lattice |(B)| ioi p = 2tt/L and 
ts = 10, and in Fig. ^ we show the angular average of ImGQ"''(t, a;;p, 0; 2^^) on the same 
lattice and again for = 10. Comparing these graphs with those for the scalar product of 
the currents, one sees that the matrix elements show here a larger variation with distance, 
and thus this should provide a stronger test of the validity of the OPE. 

In Fig. we compare the results for Re G^i^ {t, x; p, —p; 2ts) with the OPE perturbative 
predictions. Here, as always in this Section, we use the RG-improved perturbative method 
to compute 'g, using the four-loop expression for the j3 function. As we explained at length 
in the previous Section, no significant difference is observed if one uses the finite-size scaling 
method, or improved lattice perturbation theory. 

In graphs (A) and (B) we show the combination 

V{t) ^^r" [Gff(t,0;p,-p;2t,) -Gjf (-t,0;p,-p;2t,)] , (4.42) 

for two different values of p. In the scaling limit V{t) is a function of p^'^^^ and of t/^'^^'P. As 
it can be seen from the graphs our results show a very nice scaling: The data corresponding 
to the three different lattices clearly fall on a single curve. In the same graphs we also report 
the OPE prediction (^1381) 



I.e. 



f exp 

\/OP^(t) ^ 2Nu;{pr-^WnGimmt)){P,^\Jo'm -V.h). (4.43) 

Note that (p, a|c?oJo*| ~P,b) = 0, so that the corrections due to higher-order terms in the 
OPE expansion are of order t. In graph (A) and (B) we use only the one-loop Wilson 
coefficient for W^Gi^g)- There is a good agreement between the OPE prediction and the 
numerical data: quite surprisingly the agreement extends up to 2 lattice spacings. 

To better understand the discrepancies, in graphs (C) and (D) we report V{t) — V'^^^{t). 
In graphs (CI) and (C2) we use the one-loop Wilson coefficient, and in (Dl) and (D2) the 
two-loop Wilson coefficient given in Eq. ( [4.39|) . The numerical data refer to lattice |(A)| for 
(CI) and (Dl) and to lattice |(B)| for (C2) and (D2). There is clearly agreement, although 
here deviations are quantitatively somewhat large, since V{t) is strongly varying. Let us 
consider for instance the data obtained on lattice |(B)| for p = 2tt/L. If we evaluate the 



23 



matrix element of the Noether current using the numerical estimate of (t) with t = 3, 4, 5, 6 
we obtain the result with a systematic error of 2%, 14%, 32%, 27% respectively. 

In Fig. ^ we compare the angular average of Im Gq"-* (t, x; p, 0; 2^^) (cf. Fig. |) with the 
OPE prediction, by defining 

Y(r) = (4 44) 

where x*^^^'"^(r) is the OPE one-loop prediction ( [4.4(j| ). Here we use the form-factor 
prediction ([4.36|) for the matrix elements of the currents, but no significant difference would 
have been observed if the matrix elements of the currents had been determined numerically. 
Use of the form-factor prediction allows only a reduction of the statistical errors and thus 
gives the opportunity for a stronger check of the OPE. Again we observe a nice agreement 
and a very large window in which the data are well described by perturbative OPE. The 
systematic error is below the statistical one (which is approximately 10%) as soon as r > 2. 



5 Conclusions 

In this paper we have studied the OPE of the Noether currents in the two-dimensional 
(T-model. We have investigated in detail the possible problems that may arise in numerical 
applications. We can summarize the main results of this paper as follows: 

1) The relevant correlation (four-point) functions show a good scaling behaviour and 
can be determined on relatively coarse lattices (for instance on lattice |(A)|) . Lattice 
corrections of order are negligible compared to statistical and systematic errors. 
Notice that we have used the standard action, so that we could have easily improved 
the scaling behaviour by using an improved action [^] and improved operators. 

2) The running coupling constant ^(Ax) can be determined quite precisely. For instance, 
using the finite-size scaling method 'g{Ax) can be computed with an error of a few 
percent for the values of Ax of interest. This error is practically negligible in our 
applications. 

3) The main source of error is in our case the perturbative truncation of the Wilson 
coefficients. If one uses one-loop RG-improved Wilson coefficients — this also requires 
the two-loop anomalous dimensions of the operators — one can compute the matrix 
elements with a precision of 5-10%. In practice, the precision of the results is limited 
by the statistical error, which is quite large: using O(IO^) independent configura- 
tions, we obtain results with a statistical uncertainty of 10%, so that improving the 
expressions for the OPE coefficients at two loops appears of little interest. 

In conclusion the method works nicely for the cases we consider. In all cases, we observe a 
large window in which the numerical data are well described by one-loop OPE. 
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A Renormalization of composite operators 



A.l Continuum MS scheme 

In this Appendix we consider the 0{N) nonhnear a-model in d = 2 + e dimensions. In 
terms of bare fields, the action is written as 



S(a) = — [ d'^x 
9 b J 



-{d(TB{x) f - HbCTb 



(A.l) 



In order to perform a perturbative expansion, the iV-vector field ctb is parametrized as 

{TZB{x),aB{x)), 



o-b[x) = [nB[x), aB[x) ), (Tb{x) = Y 1 - 7r|j(x 

The theory is renormalized by introducing the renormalized quantities 

gB = fi'-'N^'Z.g, 

7Zb{x) = Z^^'^7t{x), 

aB(x) = Z^l^oix) = J\ - Z7z^(x), 



h 



B 



(A.2) 

(A.3) 
(A.4) 
(A.5) 

(A.6) 



where = (47r)~'^/^/r(l + e/2). Such a factor is introduced to implement naturally the MS 
prescription. The renormalization constants defined above are known to four-loop order in 
perturbation theory [Q,^,^. Here, we only need their one-loop expression: 



Z, = l 



iN-2) 
27re 



9 + Oig'), 



27re 



9 + 0{g'). 



(A.7) 



The /3-function and the anomalous dimension '^{g) of the field are obtained from the renor- 
malization constants: 



P{9) = eg 



1 + 



9 dZg 
Zg dg 



l{9) 



mdz 

Z dg 



(A.8) 



We also consider the renormalization of composite operators. In general, a renormalized 
composite operator can be written as 



[A] 



MS 



[X 



X] 



(A.9) 



B 



where the S's are unrenormalized composite operators, that is products of the renormalized 
fields tt's, a's and of their derivatives. Which operators B appear in the r.h.s. of Eq. 
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)? The naive answer is: all operators that transform under 0{N) as A and that satisfy 
dim [5] < dim [A]. This answer is wrong because of the magnetic term in Eq. ( |A.1| ) that 
breaks explicitly the 0{N) symmetry.^ The correct answer has been given in [^,|6^. If 
we define the 0{N) noninvariant operator 

a(x) = —(d^cTB + hB) = - (9V + Z,Z-^h) , (A.IO) 
as a ^ " ' 

then we must also include all products of the type ■ B' where B' has the correct 0{N) 
transformation properties and satisfies dim[i?'] + 2k < dim[yl]. Note that a{x) can be 
rewritten as 

6 S 

a{x) = hBOB - {da-Bf + Qbt^b ■ , (A.H) 

OTZb 

so that, on-shell and in the limit /i — 0, we recover an 0{N) covariant mixing matrix. In 
other words, matrix elements between physical states have the correct 0{N) transformation 
properties. 

Finally, we note that, if the MS scheme is adopted, the renormalization constants are 
dimensionless so that the operators B appearing in the r.h.s. of Eq. ( |A.9| ) satisfy dim[i?] = 
dim [A]. 

In the following we report the renormalization constants for a few cases of interest. 
A. 1.1 0{N) invariant operators of dimension 2 

We consider 0{N) invariant operators of dimension 2. There are two such operators, 
dfj^cr ■ di^cr and (dcr)'^. For the renormalization we must also consider the operator a{x) 
defined in Eq. ( [A.IOI ). They renormalize as follows:0 

[d^cT-d,a]j^ = zffd^cT-d,cT + zff{dcTf5,, + zff^a5^,, (A.12) 

[(^^)']m5 = zf,'\dcTf + zffa, (A.13) 

Nm5 = zff\dcTf + zffa. (A.14) 
A particular combination of these operators is the energy-momentum tensor 



1 „ 1 



_ = Zt,t- d,a- ■ d,cT - -6,,{day ] . (A.15) 



For this particular set of operators, the renormalization constants can be directly related 
to Z and Zg. Indeed, one first observes that T^i, is conserved. Therefore, it is finite when 
it is expressed in terms of bare fields. It follows: 

Ztt = Zff^ = ^. (A.16) 

^9 



A similar problem is encountered in nonabelian gauge theories. In order to renormalize gauge-invariant 
operators, gauge noninvariant quantities must be considered. 

^■^We adopt the following convention in reporting the renormalization constants. Given a basis 
{Qi, . . . ,Qi\ of operators of dimension d and isospin s, their renormalization constants are defined by 

[Qi\MS = 1^1 ^ij ' ''Qj- 
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Then, one notices that differentiation of the action ( |A. 1| ) with respect to g (and keeping 
d, h and tt constant) gives a finite operator, which is a combination of a and (<9cr)^. For 
g = this combination reduces to (5cr)^. Therefore, it can be identified with [{dc 



The renormalization constants immediately follow: 



7(2,0) 
^22 



7(2,0) 

^23 



1-4(1) -^«<^^'' 

-^gi-\ogZ = -^g + 0{g') 



Z/dg 



27re 



V Ja/5- 

(A.17) 
(A.18) 



A second relation is obtained by using the fact that the equations of motion do not need 
any renormalization. It follows that the combination 



is finite. Combining this result with what we obtained above, we conclude that 

Z\ 1 



'(2,0) 

32 — 


z 


7(2.0) 

~ ^22 


d 

^ dg 


'(2,0) _ 
33 — 


1 


7(2.0) 
" ^23 





2Txe 



9 + 0{g^), 



(A.19) 



(A.20) 



Z\r' = -^-Z}>r' = ^{l+g-^^\ogZj=l + ^^g + 0{g'). (A 
The remaining renormalization constants can be obtained by using Eqs. ( |A.12| - |A7TB| ): 



7(2,0) 
'12 



7(2,0) 

^13 



- z 



A2,0) 
'22 



\ d_ _ 



Aire 



9 + 0{g' 



7(2,0) 

"23 



^-9i-iogZ = -^g + 0{g') 



2Z„''dg 



A-nt 



(A.22) 
(A.23) 



Note that Eq. ( |A.11|) holds also for the renormalized operators. Indeed, using the expres- 
sions of the renormalization constants reported above, we can show that 



MS 



(x) = ha{x) - [{d(Tf]^{x)+ g-rz 



6S 
67r(x) 



(A.24) 



In the limit h and on-shell, Eq. (|A.24 ) allows us to express matrix elements of 
the 0{N) noninvariant operator [a];^ in terms of those of the 0{N) invariant operator 

Finally, it is possible to compute the trace of the renormalized energy-momentum tensor 
that gives the conformal anomaly. Using Eq. (|A.15|) we have 



[Til + T22] 



MS 



lZTT-{d(Tf. 

9 



Then, using Eqs. (|A.13|) and ( |A.14|) , we can rewrite 



{d<Tf = { ^{g) [a] 



MS 



1(9) [idar] 



MS 



(A.25) 



(A.26) 
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where we have expressed dZ/dg and dZg/dg in terms of f3{g) and 7(5') using Eqs. 
We finally obtain 



A. 1.2 Antisymmetric operators in two 0{N) indices 

We shall now consider operators that are antisymmetric 0{N) tensors with two indices. 
Obviously, there is no such operator of dimension zero. The only antisymmetric operator 
of dimension 1 is the Noether current corresponding to the 0{N) symmetry: 

fj>^^[a-dy-a%cr^]^ = Z'^''^')^ {a'^Oy - a%a^) . (A.28) 



Conservation of j'^^ implies 



Z(^'i) = ^. (A.29) 

There are three operators of dimension 2 with the same symmetry. All of them can be 
expressed in terms of dfj^j^^. Therefore, they renormalize as the current itself. 

A. 2 Lattice 

We now consider the lattice action If one neglects lattice artifacts of order a^log^a, 

lattice correlation functions of tt and a fields differ by a finite renormalization from those 
computed in any other scheme. Here we want to relate lattice correlations to the corre- 
sponding continuum ones computed in the MS scheme. The finite renormalization is given 
by 

gi = Z^g, (A.30) 
77. = (A31) 
a, = (A.32) 

where the lattice fields appear in the left-hand side, the renormalized continuum fields in 
the right-hand side. The renormalization constants Z^ and Z^ have been computed to 
three- loop order in perturbation theory ||3^,^. We shall only need the one-loop result: 



= l + ^9L\ogi^)--g, + Oigl), Z^ = 1 + ^gM^) +Oigl). 

(A.33) 

Analogously, we can relate MS -renormalized composite operators to the lattice correspond- 
ing ones. As in the continuum, besides 0(A) invariant operators, we must also consider 
operators involving a^, cf. Eq. (|3.14|) . However, also on the lattice, 0(A) invariance is re- 



covered when considering matrix elements between on-shell states. A second complication 
in the lattice approach is the loss of Lorentz invariance: as a consequence, additional opera- 
tors that are only cubic invariant must be considered. Below we report the renormalization 
constants for a few cases of interest. 
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A. 2.1 0(A^)-invariant operators of dimension 2 

On the lattice we have:[3 

+Z^t'^6,AdrT)l + Zi^i'^'^6,^a^^, (A.34) 
'^M. H<^)']ms(^) = Zi^''^5,.{d,<T)l + z!^t'^5,,(dcT)l + 

(A.35) 

[{dcTf]^{x) = Zii'''\da)l + Zii'''^a^^, (A36) 

Mm5(^) = ^4?'°^(^'^)^ + ^4?'°^«,^- (A.37) 

Note the presence of an additional operator, which is cubic but not 0{N) invariant. At 
one loop we obtain:]^ 

= l-^^.log(^) + V + 0(,i), (A38) 

Z^r^ = {l-l)9L + 0igl), (A.39) 

Z^m^Z^m ^ -^,,log(^)-^,. + 0(,i), (A.40) 

Zir-ZT' = ~H\oz(f^)-^u + 0(A), (A.41) 

4^'°' = -^S.log(^^)H-(iV-l)(±-i)<,, + Otoi),(A.43) 
= 3^SLlog(^^) + isz, + 0(ji), (A.44) 
4'''°' = l + ^9Llog(^^) + igi + 0(9i), (A.45) 

and ^'^2^^'°'' = + Z^-^"^'^^ . Finally, for the energy- momentum tensor, we have 

9l Ol 



^^In the following equations we have not reported an additional contribution proportional to where 
1 is the identity operator. This term does not contribute to connected correlation functions, and for this 
reason, we have not written it explicitly. 

Some of the renormalization constants given below have been computed in Ref. [^6[. We report them 
here to correct a few misprints. 
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with 



= l+(^l-\^9L + 0{gl), (A.47) 

= {l-il)9L + 0{gl), (A.48) 

4r = {l-\)9L + 0{9l), (A.49) 

4r = -^9L + 0igl). (A.50) 



Note that, since [T^pj-p^ is conserved, 

da 



Z^f'""^ = 0. (A.5i: 

9l 



Finally, notice that we can eliminate using the lattice equations of motion 

1 X clatt „ 

«^ = ha^ + . aV,. + -n,^^ . ^ - ^ . (A.52) 



However, — cr^. ■ d'^CTx is different from {dcr)'^ appearing in all previous formulae. The two 
operators are related by a finite renormalization: 

-(T,-d'a, = C,(d(T)l + C2ai (A.53) 



At one loop 



Then, considering only connected correlation functions, on-shell and for /i ^ 0, we have 

c^' = -T^(dcT)l. (A.55) 

Using this relation, we can compute Z^g^'^"* appearing in Eq. ( |2.10D . Explicitly 

yL,{2,0) _ yL{2,0) Cl vL(2,0) { A ^fi\ 

A. 2. 2 Antisymmetric operators in two 0{N) indices 

Let us consider operators which are 0{N) antisymmetric tensor with two indices. As we 
discussed in the continuum case, there is only one relevant operator of dimension 1, the 
0{N) current. Since the 0{N) symmetry is preserved on the lattice, there is a lattice 
discretization, given in Eq. (p.7|) , that is exactly conserved. It follows that 

j;\^)=jif, (A.57) 

i.e. Z^(i'i) = 1. 

As in the continuum there are three operators of dimension 2. If on the lattice we express 
them in terms of <9~jy'"'', then also these operators do not need any renormalization. 
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B OPE for 0{N) Noether currents 

In this Appendix we report the exphcit one-loop expressions of the OPE coefficients re- 
ported in Sec. |^. 

B.l Continuum 

At one-loop order the Wilson coefficients appearing in Eq. ( p.3|) are the following: 



N -1 



N -2 



N -1 



Attx"^ 



27C 

N 



1 - 



^(7 + log(^x)) 
- 2 



2^ ^ ( 7 + log(/xx) + - 



Wi{x) 

W2ix) 

Wsix) 
Wiix) 
W^ix) 
W,{x) 

Wjix) 
W,(x) 



N -2 



iV-2 



9 + Oig' 



2tx 
N -2 



271 



9 + Oig' 
9 + Oig' 



N — 2 

1 - -^9il + log(/^x)) + 0(/) 
3A^- 5 



9 + Oig' 



\ ^-2 N-3 . , , 

1 ■■ 9 77—9 (7 + log(^a;)) 



1 
2 

N-1 

An 
N-1 

Air 



Att 



2n 



+ 0{g'), 



g + 0{g' 



g{^ + \og{^ix))^0{g^ 



The coefficients appearing in Eq. ( |3.11| ) are given by: 

N -2 



Uoo{x) 
Ui{x) 

U2(X) 



2tt 
N-2 



g + o{g' 



An 
N-2 



g + o{g'), 



Att 



g + o{g' 



N-2 . , , N-6 ^, 2^ 
(7(7 + log(/ix)) + — g + 0{g^) 



27r 
A^-2 



An 



An 



g + o{g' 



(B.l 

(B.2 
(B.3 
(B.4 
(B.5 
(B.6 
(B.7 

(B.8 
(B.9 

(B.IO 
(B.U 
(B.12 
(B.13 
(B.IA 



B.2 Constraints on the OPE coefficients 

In this Section we want to derive the constraints due to the current conservation on the 
coefficients of the OPE. First, we need the Ward identity related to the 0{N) invariance 
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in the presence of a magnetic term h. A simple calculation gives: 



Then we need the OPE of j^\x)a' (0). The leading term for x ^ has the form 

X 



:b.i5) 



X 



(B.16) 



Using Eq. (|B.15| ) and noticing that a°'{x)a^{0) ~ 0(1) for x — >• 0, we have df{fix; g)/dx^ = 
0. Thus, /(/ix; g) is a function of g only. But the Wilson coefficient satisfies the RG equation 

Thus, if it is independent of x, and therefore of /x, it is also independent of g. A simple 
calculation at tree level gives then 



(B.18) 



The same result has been obtained in using the canonical formalism. 

We now consider the OPE of the scalar product of currents. Using the Ward identity 
( P.15|) and Eq. (|B.18| ) we have for x 



^ ab ^ ^ 



hN ~lx 



V^(0), 



g 27r X" 



(B.19) 



where we have discarded contact terms. Then, using Eq. ( |A.24| ) and discarding again 
contact terms, we obtain for x — 



2 ^ V V / 2'Kg X 

ah ^ 



This equation implies the following relations on the Wilson coefficients: 
x^d^Wo^^p{x) = 2Xf,Wo^f,p{x), 



2x' 



d 

dx"^ 
d 



[Wi{x) + W2{x) + 1^3(3;)] = 21^1 (x) - W2{x) + 2W^3(x) 



^ [^3(X) + W^{X)] = -Wi{x) - 1^3(X), 



dx"^ 

2^'^ [^5(2;) + W,{x)] 



(B.20) 

(B.21) 
(B.22) 

(B.23) 



+ Yg t'^^^) + WM - W,{x), (B.24) 



2x^^ [W^{x) + W,{x)] = ^^g + ^iig)Wsig) - W-r{x). 



(B.25) 



In the derivation we used Eq. ( [A. 27] ) for the trace of the energy-momentum tensor. 
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Now let us consider the antisymmetric case. Using the Ward identity ( [B.15 ) and Eq. 
( P.18|) , we obtain for x — 

E^/^^T(^)^'^^(0)-5MJr(^)jr(0) = (^'^V^(O) - 5^V(0)) = (0). 

(B.26) 

Again, contact terms have been discarded in the derivation. Using this relation, we obtain 
the following constraints on the Wilson coefficients^ 



d 

x^-^ [Uoo{x) + Uoi{x) + Uo2{x)] = Uoi{x) + Uo2{x), 



2x' 



_d_ 

d 



Uq2{x) = -Uoi{x) - Uo2{x), 



+ Uo2{x)] = -Uo2{x), 



X 



dx"" 



2x' 



d 



dx'^ 



[U2ix) + Uo2ix)] = -Uoi{x) - Uo2{x), 
-U2{x) + Uoo{x) + Uoi{x) + Uo2{x)] 



-2U2{x) - Uoo{x) + 2Uoi{x) + 2f/o2(x) 
iV - 2 

Uo2{x) - U2{x) = — g. 

Zn 

Using Eqs. ( p.27|) and ( p.28|) we obtain immediately 

d 



dx' 



[Uooix) + Uoiix) + 3Uo2ix)] = 0, 



(B.27) 
(B.28) 
(B.29) 
(B.30) 

(B.31) 
(B.32) 

(B.33) 



which implies that this combination is x and fi independent. By making use of the RG 
equations (|3.12| ) one proves that this combination is determined uniquely by its one-loop 
value. Then, using the results of the previous Section, we obtain: 



- 2 

Uooix) + Uoiix) + 3f/o2(x) = g. 



TX 



(B.34) 



Thus, using (p3.34| ) and ( |B.28|) , Uoo{,x) and Uoi{x) are uniquely determined by Uo2{,x). 
Moreover, using Eqs. ( [B.34| ) and ( |B.32| ) one immediately verifies that Eqs. ( |B.30| ) and 
are equivalent to Eq. {WM - 
Finally, consider ( [B.29| ). We will now show that this equation provides the two- loop 
estimate oiUo2{x). Indeed, since 



(x)] = 



(B.35) 



16 



Equations ( |B.27| )- (|B.31|) have been derived in Ref. j?^. Eq. ( |B.32|) is due to the matching of the 
terms proportional to 9^j^''(0). It is also a simple consequence of Eqs. (|B.2q ) and ( B.30 ) and of the RG 
equations. 
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and Ui{x) + Uo2{x) = 1 + 0(5'), cf. previous Section, we have 

Ui{x) + Uo2{x) = 1 - (/9olog/ix + ao)g - (/5ilog/ix + ai)g^ + 0{g^) 



[BM) 



where P{g) = —g^J2k=ol-^kg'', and Oq, cti are constants that are not fixed by the RG 
equation. Plugging this expression into ( p.29|) , we obtain 



Uo2{x) = ^g + ^g' + 0ig^) 



N-2 N-2 



-g 



g' + o{g% 



^B.37) 



Of course, the result at order g agrees with the expression reported in Sec. p.l| . Corre- 
spondingly we obtain 



UM = ^g-^g' + 0{g') 



U2(X) 



27r 
iV- 2 

N-2 
An 



g 



47r2 
N-2 

" 87r2 
N-2 



g' + oig' 



9 + ^7zr9' + 0{g'). 



(B.38) 
(B.39) 
(B.40) 



Let us finally mention that in Ref. |]70| it was argued that the functions Uoi{x) are one- loop 
exact, in the sense that there are no corrections of order g^, k > 2. As we have shown 
above and it has also been recognized by the author, ^'^ this is inconsistent with the RG 
equations. 



B.3 Lattice 

At one-loop order the Wilson coefficients appearing in Eq. ( p.l3| ) are the following: 

A^- 1 " 



W[{x) 



W,\x) 



Wiix) 



W^{x) 

W^{x) 
W}{x) 



5 



fj.p 



in 
N-1 



1 - ^V^^7l(27 + log(32x2)) - ]gL 
An A 



N-2 , , . 9^ X 1 ■ 

1 - — g (27 + log(32x2) + l) - -gL 



An 



N-2 



An 
N-2 



gi + 0{g 



2n 
N-2 



2n 



9L + 0{gl) 
9L + 0{gl) 



N-2 



An 



gL{2^ + log(32x^)) +[---]gL + 0{g 



1 1 



n 



0{gl), 
0{gl), 



+ 0{gl)iB.Al) 

(B.42) 

(B.43) 

(B.44) 

(B.45) 

(B.46) 
(B.47) 



^^In Ref. it was also shown that, even though the expressions for the Wilson coefficients were 
incorrect, one could still modify the argument so that the main result (existence of a conserved charge) of 
Ref. [[70| remains true. 
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W}{x) = 0{gl), 

Wlix) = (l-§^)9L + 0{gl), 

W,\x) = i{l-^,,-^,,(27 + log(32x^))-i,,}+0(,i), 

Wr'ix) = ^^g^ + Oigl), 

= ^^9L{27 + logi32x')-7r)+0{gl). 



B.48) 
B.49) 

B.50) 

B.51) 

B.52) 

B.53) 

B.54) 
B.55) 
B.56) 
B.57) 
B.58) 



At one-loop order the Wilson coefficients appearing in Eq. ( p.l6| ) are the following: 

N — 9 

UU^) = -^9L + 0{gl), 

T , . N -2 ^, 

U^ii^) = -^^9L + 0{gl), 

N — 2 

U^2i^) = -^9L + 0{gl), 

UHx) = l-^^^gL{2-f + \og{32x^))-^gL + ^^gL + 0{gl), 
N — 9 

U^{^) = -^9L + 0{gl). 

C Evaluation of the running coupling constant 

The determination of the running coupling constant is a key ingredient in the application 
of RG-improved perturbation theory to any asymptotically free theory. If we use the lattice 
OPE, the coupling is the bare gi, one of the input parameters of our numerical calculations. 
However, perturbation theory in gi is poorly behaved, so that one expects a poor agreement 
with the numerical data. It is known that it is much more convenient to use perturbative 
expansions in the MS scheme. Perturbative coefficients are smaller, so that truncations in 
the number of loops give smaller systematic errors. For these reasons, it is important to 
relate the MS coupling gM~s{^i) to the bare coupling g^. 

Here, we shall outline several different procedures — all of them are exact in the contin- 
uum limit — and we shall compare their efficiency. We consider the following methods: 

(A). The naive perturbative method. 

In this scheme, one computes the continuum coupling as a function of the lattice 
coupling by matching the continuum and the lattice perturbative expansion of some 
physical quantity, e.g., of the two-point correlation function. At /-loops one obtains 
a truncated series of the form fi'M5(/ia) = S'l + Yl^k=2 '^k{^cb)g\- We shall call g^^{fj,a) 
the value obtained in this way. The relevant perturbative expansions are known to 
three loops . 
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(B). The RG improved perturbative method. 

The idea of this method — and also of those that will be presented below — is to com- 
pute gjisifJ') starting from some quantity that can be computed numerically at the 
given value of the bare lattice coupling constant. 

In this case, we consider the RG prediction for the mass gap m in the MS scheme: 



(C). 



where 



A- 



A/5 



fx e 



exp 



11 /3i 



dz 



(C.l) 



(C.2) 



and I3{g) = —g'^^PuQ^ is the MS /3-function. The constant Cn is not known for a 
general theory. For the two-dimensional cx-model it has been computed P9| , [50[] using 
the thermodynamic Bethe ansatz. The result, in the MS scheme, reads: 



C 



N 



r(i + (iv-2) 



(C.3) 



The method works as follows: For a given value of the lattice coupling ql, compute 
numerically (for instance, by means of a Monte Carlo simulation) the mass gap ma. 
Then, fix fia and solve numerically Eq. ( |C.1[ ), obtaining gjYg{fi/m). Note that, since 
the /^-function is known only to a finite order in perturbation theory, we have to 
substitute the function Ajg^(/i, g) with its truncated perturbative expansion. There 
is some arbitrariness in this truncation. We shall make the simplest choice 

1-2 



1 + 5^ A,/ 



k=l 



(C.4) 



where the coefficients are obtained by expanding perturbatively Eq. ( |C.2| ). Equa- 
tion ( |C.4| ) gives the /-loop approximation of the A-parameter. The solution of the 
corresponding Eq. ( |C.1| ) will be denoted as g^^''{lJ'/'m)- The perturbative expansion 
of the MS /3-function is known to four loops p4|. 



The finite-size nonperturbative method. 

due to Liischer [^], was initially tested in the two-dimensional 0(3) 



This method 
(T-model HUT 



Recently, it has been successfully employed in the computation of the 



A-parameter in quenched QCD [|T^. The idea is to consider the theory in a finite 
box and to define a "finite-size scheme" in which the renormalization scale is the 
size of the box. For the a-model, Ref. introduces a coupling gn^l/L) defined as 
follows:|3 

MUL) = (C.5) 



^^This definition is by no means unique. For instance, one could also use gji{l/L) = [ni{L)L]'^ / {N — 1), 
where m(L) is the inverse of the second- moment correlation length on a square lattice of size L/a. The 
corresponding universal finite-size scaling function — i.e. the function that gives the correspondence between 
(7fl(l/L) and mL — has been determined numerically in |4q|. 



36 



where m{L) is the mass gap in a strip of width L. Standard finite-size scahng theory 
indicates that gji is a universal function of mL, where m is the infinite-volume mass 
gap. Such a function can be computed nonperturbatively by means of Monte Carlo 
simulations with a good control of the systematic errors. If we set fi = 1/L, gji defines 
a running coupling constant that is a function of fi/m. The function gji{l/L) can also 
be computed in perturbation theory. This provides the connection between and 
any other perturbative scheme. This immediately defines the Z-loop approximation 
to the MS coupling: g^^(fi) = gR^fi) + X]fc!=2 '^fcfl'M/^)- The perturbative expansion 
of gji{l/L) is known to three-loop order fT^ . 

(D). The hybrid method. 

The finite-size scaling method does not provide — at least in the implementation of 



Ref. ||6T[] — the coupling gjjsifJ,) for any /i, but only on a properly chosen mesh of 
values, say {fii}i=i,.... The next step consists in interpolating between these values 
using RG equations. Since the mass gap is a RG-invariant quantity, at order /, we 
may require 

^^^{^^^,ffMs{^^^)) = Ag^(/x,^7^(/.)). (C.6) 

Using gjjs(f^i)^ '^^^ '^^^ then obtain gjjg{jj) for any given ^. We shall call g^^{^) 
the running coupling obtained by this procedure. 

(E). The improved- coupling method. 

Method |(A)| does not work well because lattice perturbation theory is not "well be- 
haved": Perturbative coefficients are large, giving rise to large truncation errors. 
Parisi |62,^| noticed that much smaller coefficients are obtained if one expands in 



terms of "improved" (or "boosted" ) couplings defined using "short-distance" observ- 
ables. In the cr-model one can define a new coupling in terms of the energy density 

4 

which is then related to gjjsil^) perturbatively. At order I, we can write g^^{afi) = 

gE+Yl^k^2 '^k{l^^)9E- practice the method works as follows: Given gi, one computes 
numerically gE', then, one uses the previous perturbative expansion to determine the 
MS coupling constant. This method is expected to be better than the naive one. 
Indeed, one expects |cf(/ia)| <C |cfc(/ia)|, so that truncation errors should be less 
important. The perturbative coefficients cf can be computed up to / = 3 using the 
results of 



Notice that the list above is by no means exhaustive. For instance, an alternative nonper- 
turbative coupling may be defined using off-shell correlation functions: 

9Ri^^ = n = e^ii^i^Ai^^ . (C.8) 
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m 


j_ 




9jTS ^^'> 






0.5372 


0.00071(11) 


0.0097(15) 


0.5870 [-3] 


0.5892[-l] 


0.574[-5] 


0.5902[-28] 


0.5747 


0.00143(11) 


0.0195(15) 


0.6321 [-4] 


0.6351 [-1] 


0.617[-6] 


0.6354[-31] 


0.6060 


0.00237(15) 


0.0323(20) 


0.6703 [-5] 


0.6736 [-2] 


0.652[-6] 


0.6737[-16] 


0.6553 


0.00478(15) 


0.0652(20) 


0.7312[-7] 


0.7362 [-2] 


0.708 [-8] 


0.7364[+5] 


0.6970 


0.00794(19) 


0.1083(25) 


0.7835 [-9] 


0.7900 [-3] 


0.755 [-12] 


0.7903[+6] 


0.7383 


0.01231(15) 


0.1678(20) 


0.8361[-11] 


0.8437[-4] 


0.800[-16] 


0.8436[-13] 


0.7646 


0.01589(15) 


0.2166(20) 


0.8701 [-13] 


0.8788 [-5] 


0.828 [-20] 


0.8777[-35] 


0.8166 


0.02481(22) 


0.3382(30) 


0.9382[-16] 


0.9486 [-7] 


0.881 [-28] 


0.9434[-99] 


0.9176 


0.04958(38) 


0.6759(51) 


1.0742 [-26] 


1.0852[-13] 


0.975 [-45] 


1.0623 [-275] 


1.0595 


0.1033(6) 


1.4082(77) 


1.2743 [-47] 


1.2895[-27] 


1.090 [-73] 


1.2135[-593] 


1.2680 


0.2092(5) 


2.8519(67) 


1.5886[-96] 


1.59626[-680] 


1.217[-108] 


1.3863[-1056] 



Table 4: The MS running coupling constant. We use here several different methods as 
explained in the text. The errors on the second and third columns are statistical. 



Something similar has been proposed in Refs. [|73|-[75|], with the purpose of computing the 
QCD A-parameter. This approach opens the Pandora box of possible definitions of the 
running coupling in substitution of Eq. ( U.8 ). A scheme that has been intensively studied 
in the context of QCD employs the three-gluon vertex (see Refs. jT^, [76H8T[] ) . 

Let us compare the different methods. In Tab. ^ we compare the procedures |(A)| , |(B) , 
(C)l , and |(E)| . In the first column we report a collection of values of qr. A subset of the 
values given in the table have been considered for the first time in Ref. [^. Later, the 
mesh was enlarged by Hasenbusch |p2|]. For these values of qr, Hasenbusch computed the 



corresponding value of mL = m/ jj, which is reported in the second column. Note a pecu- 
liarity of the finite-size approach: usually, one fixes /i/m and then determines the running 
coupling constant. Here, the running coupling constant is fixed at the beginning and the 
value of the scale is determined numerically. In the third column we report the scale in 
lattice units for gi = 1/1.54, the value of the lattice coupling at which we have done most of 
our simulations. The results are obtained by using [ma)~^ = 13.632(6). The error reported 
there corresponds to the error on m/n, the error on (ma) being negligible. In column 4 we 
report the estimate of g^^(jj?j obtained by using and three-loop perturbation theory [^ 



In brackets we report the difference g^^{fi) — g^^{fi). In the next column we report the 



fs,3 , 



'MS 



four-loop coupling g^-^^[ii) obtained by using the value of m/fi given in the second column. 



Again, in brackets we report g^^^ijJ') — 5'il^^(/^)- In the last two column we report the 

I, 1-3— r 

results obtained by using three-loop lattice perturbation theory ||3^, . In the fifth 
column we use gi = 1/1.54 as the expansion parameter. In the sixth column the improved 
coupling defined by Eq. ( |C7|) is used. The connection with the bare coupling is obtained 
by using the perturbative expressions given in Ref. 
has been evaluated in a Monte Carlo simulation at gi 
statistics A^'stat = 10000, yielding gs = 0.768133(49). 
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The relevant expectation value 
1/1.54 on a lattice 128 x 256 with 
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9ms 


3ms ^^'> 


hybr,i(/ \ 

SjTS ^^') 


1 


1.18422 


1.171(2) 


1.1804(6) 


2 


1.42282 


1.403(3) 


1.417(1) 


3 


1.62532 


1.598(4) 


1.617(1) 


4 


1.81923 


1.783(5) 


1.809(2) 


5 


2.01713 


1.970(7) 


2.003(2) 


6 


2.22905 


2.167(9) 


2.211(3) 


7 


2.46722 


2.385(12) 


2.443(4) 


8 


2.75208 


2.637(16) 


2.717(6) 



Table 5: The MS running coupling constant. In the first column we report the RG coupling 
obtained using (ma)"^ = 13.632(6) at l/gi = 1-54. The coupling reported in the third and 
fourth column are obtain using the hybrid scheme: the columns differ in the choice of the 
"integration constant", see text. The reported error is due to the error on /ij appearing in 
the l.h.s. of Eq. (|C6|), see Tab. 0, second column. 



In Tab. ^ we compare, on a broad range of scales, the outcome of RG-improved per- 
turbation theory and the interpolation procedure which we called "hybrid." In both cases 
four-loop perturbation theory is used. The two hybrid couplings referred to as A and B 
differ in the "boundary condition" for the RG interpolation, that is in the value used in 
the left hand side of Eq. (|C.6|) . We used respectively, gj^ = 1.074254 for m//i = 0.04958 
and gjfg = 1.588556 for m//i = 0.2092, obtained using the finite-size scaling method, cf. 
Tab. ^. In all cases we fixed {ma)~^ = 13.632(6). 

What do we learn from this comparison? First of all, lattice (naive) perturbation theory 
(sixth column of Tab. ^ is a very bad tool. Even at energies as high as 50 times the mass 
gap gjYsifJ') is affected by a ~ 5% systematic error. However, it is reassuring that the 
expansion tells us its own unreliability. Indeed, the observed discrepancy is of the order 
(at most twice as large) of the difference between the two- loop and the three- loop result. 
The perturbative expansion in terms of the improved coupling is much better. The results 
are quite precise up to ~ 10m. For smaller values of the discrepancy increases, but it 
is nice that it is again of the order of the difference between the two- and the three-loop 
result. Perturbative RG supplemented with the prediction ( |C.3| ) gives results which are 
in agreement with the nonperturbative ones obtained using the finite-size scaling method 
within a few percent for all the energy scales given in Tab. |[ The accuracy remains good 
(if the comparison is made with the "hybrid" procedure) also for scales of the order of the 
mass gap. 

Up to now we have discussed the MS scheme and how to obtain the value of the MS 
coupling. However, as we already mentioned above, reasonably good results can also be 
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obtained if we use the coupling gE- In this scheme we introduce the parameter 



-Pi /Pi 



exp 



1 



l3ogE 



dt 



1 



1 



/3i_ 

j3E{t) ^ Pof" I3lt 



(C.9) 



where Pe{9e) = -9lY.k=oPk 9%- 9l = f{gE), then Pe{9e) = P'^'\f{9E))/ f'{gE)- The 
mass gap is invariant and thus m = Cj^^e^e^cLiQe)) where 



C 



N,E 



l/(7V-2) 



r(i + ]vb) 



25/2 



exp 



TT 



_4(iV-2)_ 



(C.IO) 



References 

[1] M. Bochicchio, L. Maiani, G. Martinelh, G. C. Rossi, and M. Testa, NucL Phys. B262, 
331 (1985). 

[2] G. Curci and G. Veneziano, NucL Phys. B292, 555 (1987). 

[3] S. Caracciolo, G. Curci, P. Menotti, and A. Pehssetto, NucL Phys. B309, 612 (1988). 

[4] S. Caracciolo, P. Menotti, and A. Pelissetto, NucL Phys. B375, 195 (1992). 

[5] G. Martinelli, C. Pittori, C. T. Sachrajda, M. Testa, and A. Vladikas, NucL Phys. 
B445, 81 (1995), |E^at/94110iq . 

[6] M. Gockeler, R. Horsley, H. Oehich, H. Perlt, D. Petters, P. E. L. Rakow, A. Schafer, 
G. Schierholz, and A. Schiller, NucL Phys. B544, 699 (1999), |hep-lat/ 9807044 . 

[7] K. Symanzik, NucL Phys. B190, 1 (1981). 

[8] M. Liischer, R. Narayanan, P. Weisz, and U. Wolff, NucL Phys. B384, 168 (1992), 
|hep-lat/ 92070091 . 

[9] M. Liischer, R. Sommer, U. Wolff, and P. Weisz, NucL Phys. B389, 247 (1993), 
|hep-lat/9207010 . 

[10] M. Liischer, R. Sommer, P. Weisz, and U. Wolff, NucL Phys. B413, 481 (1994), 
|hep-lat/ 9309005 . 

[11] K. Jansen, C. Liu, M. Liischer, H. Simma, S. Sint, R. Sommer, P. Weisz, and U. Wolff, 
Phys. Lett. B372, 275 (1996), |hep-lat/ 95 12009 . 

[12] S. Capitani, M. Liischer, R. Sommer, and H. Wittig, ALPHA Collaboration, NucL 
Phys. B544, 669 (1999), |hep-lat/9810063| . 

[13] M. Liischer, Advanced lattice QCD, in Les Houches Summer School "Probing the 
Standard Model of Particle Interactions", edited by F. David and R. Gupta, 1998, 
hep-lat/ 98020291 . 



40 



[14] M. Liischer, S. Sint, R. Sommer, and H. Wittig, Nucl. Phys. B491, 344 (1997), 
|hep-lat/9611015 . 



[15] C. Dawson, G. Martinelli, G. C. Rossi, C. T. Sachrajda, S. Sharpe, M. Talevi, and 
M. Testa, Nucl. Phys. B514, 313 (1998), |h^at/9707009| . 



[16] G. C. Rossi, Talk given at Cortona 1998, Cortona, Italy, unpublished, |hep-lat/9811009 



[17] M. Testa, Nucl. Phys. B (Proc. Suppl.) 63, 38 (1998), |hep-lat/9709044 



[18] G. MartineUi, Nucl. Phys. B (Proc. Suppl.) 73, 58 (1999), |hep-lat/98100T3 
[19] K. G. Wilson, Phys. Rev. 179, 1499 (1969). 

[20] W. Zimmerman, in Lectures on Elementary Particles and Quantum Field Theory, 
1970 Brandeis University Summer Institute in Theoretical Physics, edited by S. Deser, 
H. Pendleton, and M. Grisaru, MIT Press, Cambridge, 1970. 

[21] P. A. M. Dirac, Proc. Cambridge Phil. Soc. 30, 150 (1934). 

[22] M. D'Elia, A. Di Giacomo, and E. Meggiolaro, Nucl. Phys. B (Proc. Suppl.) 73, 515 
(1999), |hep-lat/98110T7l 



[23] M. D'Eha, A. Di Giacomo, and E. Meggiolaro, Phys. Rev. D59, 054503 (1999), [Eep^ 
I lat/ 98090551 . 



[24] M. D'Elia, A. Di Giacomo, and E. Meggiolaro, Nucl. Phys. B (Proc. Suppl.) 63, 465 
(1998), |hep-lat/ 97090311 



[25] M. D'Elia, A. Di Giacomo, and E. Meggiolaro, Phys. Lett. B408, 315 (1997), [hep^ 



lat/9705032 



[26] M. Campostrini, A. Di Giacomo, and G. Mussardo, Z. Phys. C25, 173 (1984). 

[27] S. Capitani, M. Gockeler, R. Horsley, H. Oelrich, D. Petters, P. Rakow, and G. Schier- 
holz, Nucl. Phys. B (Proc. Suppl.) 73, 288 (1999), |hep-lat/9809T7l . 



[28] S. Capitani, M. Gockeler, R. Horsley, D. Petters, D. Pleiter, P. Rakow, and G. Schier- 
holz, Nucl. Phys. B (Proc. Suppl.) 79, 173 (1999), |hep-ph/9906320 . 



[29] A. M. Polyakov, Phys. Lett. B59, 79 (1975). 

[30] E. Brezin and J. Zinn- Justin, Phys. Rev. B14, 3110 (1976). 

[31] W. A. Bardeen, B. W. Lee, and R. E. Shrock, Phys. Rev. D14, 985 (1976). 

[32] E. Brezin, J. Zinn- Justin, and J. C. Le Guillou, Phys. Rev. D14, 2615 (1976). 

[33] S. Hikami and E. Brezin, J. Phys. All, 1141 (1978). 

[34] W. Bernreuther and F. J. Wegner, Phys. Rev. Lett. 57, 1383 (1986). 



41 



[35] M. Falcioni and A. Treves, Nucl. Phys. B265, 671 (1986). 
[36] S. Caracciolo and A. Pelissetto, Nucl. Phys. B420, 141 (1994), |hep-lat/94010T5 



[37] S. Caracciolo and A. Pelissetto, Nucl. Phys. B455, 619 (1995), |hep-lat/95100T5 



[38] B. Alles, A. Buonanno, and G. Cella, Nucl. Phys. B500, 513 (1997), |E^at/9701001 
[39] D.-S. Shin, Nucl. Phys. B546, 669 (1999), [Fep^at/9810025| . 



[40] B. Alles and M. Pepe, Nucl. Phys. B563, 213 (1999), |h^lat/9906012 . 



[41] B. Alles, S. Caracciolo, A. Pelissetto, and M. Pepe, Nucl. Phys. B562, 581 (1999), 
|hep-lat/9906014| . 



[42] B. Alles and M. Pepe, Nucl. Phys. B (Proc. Suppl.) 83-84, 709 (2000), [Eep^ 
I lat/9909aTP|. 



[43] U. Wolff, Phys. Lett. B248, 335 (1990). 
[44] U. Wolff, Nucl. Phys. B334, 581 (1990). 

[45] R. G. Edwards, S. J. Ferreira, J. Goodman, and A. D. Sokal, Nucl. Phys. B380, 621 
(1992), |liep-lat/9112002 , 



[46] S. Caracciolo, R. G. Edwards, A. Pelissetto, and A. D. Sokal, Phys. Rev. Lett. 75, 
1891 (1995), [Ei^at/9411009 . 



[47] S. Caracciolo, R. G. Edwards, T. Mendes, A. Pelissetto, and A. D. Sokal, Nucl. Phys. 
B (Proc. Suppl.) 47, 763 (1996), ^ep-lat/ 950903^ . 

[48] B. Alles, G. Cella, M. Dilaver, and Y. Giindiig, Phys. Rev. D59, 067703 (1999), 
|hep-lat/9808003| . 



[49] P. Hasenfratz, M. Maggiore, and F. Niedermayer, Phys. Lett. B245, 522 (1990). 
[50] P. Hasenfratz and F. Niedermayer, Phys. Lett. B245, 529 (1990). 
[51] U. Wolff, Phys. Rev. Lett. 62, 361 (1989). 
[52] U. Wolff, Nucl. Phys. B322, 759 (1989). 

[53] S. Caracciolo, R. G. Edwards, A. Pelissetto, and A. D. Sokal, Nucl. Phys. B403, 475 
(1993), PTep-lat/9205005 , 



[54] S. Caracciolo, A. Montanari, and A. Pelissetto, Nucl. Phys. B (Proc. Suppl.) 73, 273 
(1999), |hep-lat/9809T00l 



[55] S. Caracciolo, A. Montanari, and A. Pelissetto, Nucl. Phys. B (Proc. Suppl.) 83-84, 
875 (2000), |hep-lat/9909T05| . 



[56] A. Buonanno, G. Cella, and G. Curci, Phys. Rev. D51, 4494 (1995), |hep-lat/9506TO 

42 



[57J 


M 


[do\ 


T 

J . 




iVi 


[oUJ 


i\/r 
iVi 


[bij 


iVi 


[D/J 




[63] 


G. 


[64] 


G. 


[65] 


M 


[66] 


K. 



et al., Springer, Berlin, 1982; K. Symanzik, Nucl. Phys. B226, 187, 205 (1983); P. 
Hasenfratz and F. Niedermayer, Nucl. Phys. B414, 785 (1994), [hep-lat/9308Un^ ; S. 
Caracciolo, A. Montanari, and A. Pelissetto, Nucl. Phys. B556, 295 (1999), |iep^ 
Iat/9812014| . 



[67; 

[68 
[69 
[70 
[71 
[72 
[73 

[74; 

[75 

[76 

[77; 



[7^ 



S. Hikami, Phys. Lett. 98B, 208 (1981). 
F. Wegner, Nucl. Phys. B316, 663 (1989). 

E. Brezin, J. Zinn- Justin, and J. C. Le Guillou, Phys. Rev. B14, 4976 (1976). 
M. Liischer, Nucl. Phys. B135, 1 (1978). 
M. Liischer, unpublished. 

D.-S. Shin, Nucl. Phys. B496, 408 (1997), [hep^at/9611006 . 



D. Becirevic, P. Boucaud, J. P. Leroy, J. Micheli, O. Pene, J. Rodrfguez-Quintero, and 
C. Roiesnel, Phys. Rev. D60, 094509 (1999), [Ei^h/9903364| . 



D. Becirevic, P. Boucaud, J. P. Leroy, J. Micheli, O. Pene, J. Rodrfguez-Quintero, and 
G. Roiesnel, Nucl. Phys. (Proc. Suppl.) 83-84 (2000) 159, |hep-lat/9908056| . 

D. Becirevic, P. Boucaud, J. P. Leroy, J. Micheli, O. Pene, J. Rodrfguez-Quintero, and 
G. Roiesnel, Phys. Rev. D61 (2000) 114508, |hep-ph/99l0"^ . 

D. Becirevic, V. Lubicz, G. Martinelli, and M. Testa, Nucl. Phys. B (Proc. Suppl.) 
83-84 (2000) 863, |hep-lat/9909039| . 



D. Henty, G. Parrinello, and G. Pittori, Proc. of the EPS/HEP95 Gonference, [hep- 
lat/ 95 10045] . 



G. Parrinello, Nucl. Phys. B (Proc. Suppl.) 34, 510 (1994), |hep-lat/9311065 . 



43 



[79] C. Parrinello, Phys. Rev. D50, 4247 (1994), |EiFlat/9405024 . 

[80] C. Parrinello, D. G. Richards, B. Alles, H. Panagopoulos, and C. Pittori, UKQCD 
Collaboration, Nucl. Phys. B (Proc. Suppl.) 63, 245 (1998), |hep-lat/9710053i 

[81] B. Alles, H. Panagopoulos, C. Parrinello, C. Pittori, and D. G. Richards, Nucl. Phys. 
B502, 325 (1997). 

[82] M. Hasenbusch, unpublished, cited in [jT2|. 



44 



U 

O 




10 15 

iterations 



Figure 1: Logarithm of the autocorrelation function log A{i) for the three lattices employed. 
Continuous lines refer to lattice |(A)| , dotted lines to lattice |(B)| and dashed lines to lattice 
(C) . The upper curves refer to the short-distance observable Od, d = 1. The lower curves 



refer to the long-distance observable Od for d ~ ^'^^p- d = 8 for lattice |(A)| , d = 15 for 
lattice |(B)| and d = 30 for lattice |(C) . 
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Figure 2: Estimates of lmG^^\t, x;p,p;2ts) averaged over rotations on lattice KA)| (upper 
graph, here tg = 10) and |(B)| (lower graph, ts = 6). Circles, filled squares, and triangles 
correspond to p = 27i/L, 47r/L, and Qtt/L respectively. 
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Figure 3: The scalar product of two Noether currents compared with the OPE prediction: 
graphs of R{r), cf. Eq. ( [4.33| ), obtained using MS RG- improved perturbation theory. 
Circles, filled squares, and triangles correspond top= 27c/L, 47r/L, and Gtc/L respectively. 
The data are for lattice [B], ^''^p = 13.636(10). 
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Figure 4: The scalar product of two Noether currents compared with the OPE prediction: 



graphs of Eq. i?^'***(r), cf. (f4.34 ), obtained using RG-improved perturbation theory in 



the couphng gi and in the improved couphng qe- Circles, filled squares, and triangles 
correspond to p = 2tt/L, 4:71/ L, and 6tt/L respectively. The data are for lattice 
^exp = 13.636(10). 
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Figure 5: The scalar product of two Noether currents compared with the OPE prediction: 
graphs of R{r), cf. Eq. ( [4.33| ), obtained using MS RG- improved perturbation theory. 
Circles, filled squares, and triangles correspond to p = 27r/L, 47r/L, and Gtt/L respectively. 
The data are for lattice |(A)| , ^"""p = 6.878(3). Notice the change of vertical scale compared 
to Figs. I, |. 
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Figure 6: Estimates of ImG'-^^t, x; p, 0; 20) averaged over rotations on lattice |(B)| . Circles, 
filled squares, and triangles correspond to p = 27i/L, 47r/L, and 6n/L respectively. 
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Figure 7: The scalar product of two Noether currents compared with the OPE prediction: 
graphs of S{r), cf. Eq. (|4.35| ), obtained using MS RG-improved perturbation theory. 
Circles, filled squares, and triangles correspond to p = 2n/L, 47r/L, and Gii/L respectively. 

[B]|, ^^^^P = 13.636(10). 



The data are for lattice 
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Figure 9: Angular average of IrnGjif (t,a;;p,0;20) on lattice [[B]]. Circles, filled squares, 
and triangles correspond to p = 27c/L, 47r/L, and 6tt/L respectively. 
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Figure 10: Antisymmetric product of two Noether currents for x = and t ^ 0: in graphs 
(A) and (B) we report estimates of V{t), cf. Eq. ( [4 .421 ), and of the OPE prediction (dashed 
hne) cf. Eq. ( [4.43|) . The numerical data correspond to p = 27t/L (graph (A)) 

and p = Att/L (graph (B)). Stars, diamonds, and circles refer to lattices |(A)| , |(B)| , |(C)| 
respectively. In graphs (C) and (D) we show V{t) — V'~'^^{t). Empty circles and filled 
squares refer to p = 2t: / L and p = An/L respectively. (CI) and (Dl) refer to lattice |(A) 
(C2) and (D2) to lattice [B]. 
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Figure 11: The antisymmetric product of two Noether currents compared with the OPE 
prediction: graphs of Y{r), cf. Eq. ( [4 .441 ), obtained using MS RG-improved perturbation 
theory. Circles, filled squares, and triangles correspond to p = 27i/L, 47r/L, and Gn/L 
respectively. The data are for lattice |(B)| , C,'^^^ = 13.636(10). 
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